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GENERAL DESCRIPTION OF WORK, ETC. 

The work described concerns the last period of electrochemical 
research carried out for N.A.S.A. in the Electrochemistry Laboratory. 

To gain perspective, it is stated that, over the ten years of continuous 
support of N.A.S.A. to the Electrochemistry Laboratory, the basis 
governing the direction of the work was N.A.S.A.'s interest in electro- 
chemical energy conversion, i.e., in fuel cells. For this reason, much 
of the work published earlier in the decade was with respect to topics 
such as the potential of zero charge, electrocatalysis, and porous 
electrode theory. From 1968, the emphasis in the work changed, as a 
consequence of the technological success of N.A.S.A.'s fuel cell efforts. 
The new emphasis went onto electrochemical energy storage devices. 
However, in keeping with all the work supported by N.A.S.A. in the 
Electrochemistry Laboratory, work done was entirely fundamental in 


character. 
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WORK OF DR. B. BUCK 

ADSORPTION OF DENDRITE INHIBITORS ON ZINC 
INTRODUCTION 

The formation of dendritic zinc is a major obstacle 
hindering the attainment of optimum efficiency in alkaline zinc 
batteries. Previous work 1 in this laboratory has showed the 
efficacy of quaternary amines in inhibiting the growth of zinc 
dendrites during the charging cycle in alkaline zincate solutions. 
Dendritic growth occurs during the charging cycle at cathodic over- 
potentials greater than 70 millivolts. 

In order to most effectively exploit the observation that 
some quaternary amines inhibit zinc dendrite growth, it is necessary 
to understand the parameters governing the adsorption of various 
quaternary amines on a solid zinc electrode. These parameters are, 
in general, obtainable from electrocapillary data. Unfortunately 
however, zinc is not an ideally polarizable electrode; consequently 
electrocapillary data is or even pseudo-electrocapillary data is 
experimentally unobtainable. 

Since the thermodynamic quantities (surface excesses, etc.) 
are not available to us, the only alternative is to measure the total 
amount of the quaternary amines adsorbed in the interphase as a function 
of its bulk concentration and as a function of the zinc electrode poten- 
tial. There are several experimental techniques that are possibly 
useful for such measurements. 



Two possible electrochemical methods of adsorption study 
are differential capacitance measurements and oxidation of the adsorbed 
inhibitor. The difficulty with the latter technique is that it 
necessitates an assignment of a fraction of a total current to the 
oxidative process of the adsorbate; such assignment in the case of 
a zinc electrode which has a large exchange current density is tenuous 
at best. The situation with any differential capacitance technique 
is even more arbitrary in that a model must be invoked in order to 
arrive at an amount adsorbed. The overriding difficulty of both 
methods, however, is that hydrogen is adsorbed in the potential 
region of interest. 

The situation with non-electrochemical techniques is 
marginally better. The obvious methods here are of two kinds: 
analytical determination of the change in bulk concentration of 
adsorbate and measurement of the amount of adsorbate on (or very near) 
the electrode. The difficulty inherent in the first method is the 
very low (ca. 10 ^ M) bulk concentration of the adsorbate; few tech- 
niques are available for determining quaternary amines at this 
concentration level. Two techniques are available for measurement 
of the total amount adsorbed. They are ellipsometry and radiotracers . 

PRELIMINARY INVESTIGATIONS 

Preliminary investigations with the ellipsometric method 
showed it to be unsuitable for the present problem. The essential 
difficulty was the impossibility of maintaining a specularly reflective 
zinc electrode surface. The experimental conditions were such that 
roughening of the zinc surface was unavoidable. 



The above considerations effectively limit the suitable 

methods of study to radioactive techniques of vjhich there are two 

general types. Those in which count rates are determined in situ 

and those in which the adsorbent is removed from the bulk adsorbate 

and counted externally. The latter type of technique was chosen as 

the type most amenable to the present problem. 

The technique selected was developed in this laboratory 

in the early 1960's by Swinkels, et al 2 . In essence, the method 

14 

consists of adsorbing C labeled organic compounds on a continuous 

metal tape, which is connected as a working electrode in a usual 

3-electrode potentiostat circuit. The tape electrode, when slowly 

withdrawn from the cell at a uniform speed through a slit, has only 

a quite thin (ca. 1 micron) film of solution on it. The amount of 
14 

C -labeled compound adsorbed on the tape (together with that in the 
adhering solution film for which a correction can be made) is determined 
by means of two parallel - one above and one below the tape - gas flow 
proportional counters. 

EXPERIMENTAL DIFFICULTIES 

The two experimental difficulties immediately encountered 
were insufficient sensitivity of the radiation detection system and 
grossly irreproducible solution film thicknesses. Radiation detection 
system sensitivity is improvable to a large extent by changing the 
detectors and/or their geometry. The more serious problem is that of 
film uniformity. 
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The obvious cause of the latter difficulty was, of 

course, the slit system. The slits were first fabricated of Teflon. 

The Teflon has both a low hardness and a low surface free energy. 

The low hardness and lack of dimensional stability of the Teflon is 

the cause of excessive wear when used in conjunction with the metal 

tape electrode. More serious still is the low surface free energy 

of the Teflon, the result of which is that it is not water wet whereas 
the tape electrode is; the mechanical instability thus created causes 

the film thickness to be unduly sensitive to slight pressure changes 

at the exit slit. 


EXPERIMENTAL DETAILS 

A new cell with glass exit slits was designed and built (see 
previous reports). Even with the most carefully made slits, however, 
it was not possible to prevent leakage of the solution at the slit, tape, 
solution juncture when even very low gas pressures were applied to the 
drying compartment. It is, of course, necessary to do this in order 
to sparge the solution. The ultimate solution to this problem was the 
careful and accurate regulation of a vacuum system (water aspirator 
as a source) connected to both the cleaning and adsorption cells. 

The sparging gas. He, was admitted through the drying compartment. 

Very fine metering valves (Nupro) were found suitable for control of 
the pressure and vacuum systems. 

The insensitivity of the radiation detection system was 
traced to the anode in the gas flow detectors. New anodes were 

fabricated in the laboratory from 1 mil platinum wire which was soldered 
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onto a stainless steel tube that was made part of an HN connector. 

With the new anodes and using PR-10 gas, the detectors exhibited Geiger 

regions of over 300 volts with a slope of about 2% per 100 volts. 

Since the E for the beta decay of carbon 14 is about 

150 keV, the thinnest possible windows are required for the detectors 

in order to gain sensitivity sufficient to make meaningful measurements 

on the very small amounts of radiation to be measured (of the order of 
-9 2 

10 moles). Micro (85 yg/cm ) gold foil windows are the best 
available choice of windows in as much as they are the least 
adsorptive ones obtainable. 

14 

The tape detector system was calibrated with labeled Na 2 C- 0^. 

* 

A 0.001 M solution of the Na^O^ was prepared by determinative weighing. 

A small quantity (typically 0.1 ml) of this solution was placed on a 
clean tape with a micropipette and allowed to spread. The tape was located 
between the detectors in the same position it would occupy when adsorp- 
tion were made. The apparent counter efficiency obtained was 190%. 

Note that this figure represents the total number of counts by two 

2 

detectors divided by the calculated number of disintegrations per cm 

of tape. The actual area of tape counted is indeterminate, but in 
2 

excess of 10 cm . 

In the actual adsorption studies the procedure was to clean 
the tape in the cleaning cell which contained 1 M NaOH in conductivity 
water by alternate oxidation then reduction of the electrode. The 
bright electrode was then moved through the drying compartment which 
contained dry He to the adsorption cell. In the adsorption cell the 
Zn tape electrode was potentiostated at a pre-selected potential using 
a saturated calomel electrode as reference. The tape potential was 



- 7 - 


maintained for a time sufficient to allow the establishment of adsorption 
equilibrium. The amount of time required - on the order of 15 minutes - 
was determined experimentally by plotting observed count rate vs time of 
adsorption. After the suitable period of time had elapsed, the tape, 
still under potential control, was moved from the adsorption cell to its 
position between the two gas flow detectors and counted. The count 
time was adjusted so that at least a minimum of 1,000 counts was 
recorded. 


EXPERIMENTAL MEASUREMENTS CARRIED OUT 
Adsorption measurements were made on zinc tape of 99.999% 
purity which was furnished by Cominco American, Spokane, Washington, 
in the form of a strip 2 mils thick by 500 mils wide. The adsorption 
solution was prepared by diluting a stock solution of the appropriate 

quaternary amine with 1 M K0H to give the desired concentrations. 

— 6 -6 

The concentrations: selected for study were: 10 M, 2.5 x 10 M, 

5 x 10~ 6 M, 7.5 x 10 6 M, 10~ 5 M, and 2.5 x 10~ 5 M. The cetyl trimethyl 
ammonium bromide used was furnished by Amersham/Searle with a specific 
activity of 6 mCi/mM. 

RESULTS 

Results of these experiments indicate that cetyl trimethyl 
ammonium bromide is strongly adsorbed at potentials within 200 mv of the 
open circuit potential (ca. -1.5 V vs SCE) . The amount adsorbed in this 
potential region is only weakly concentration dependent. At potentials 
more than 200 mv cathodic to the open circuit potentials the quaternary 
amine desorbs strongly. 
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The at first surprising result that an organic cation 
is desorbed by a negatively charged electrode is at second glance 
quite reasonable. At the large potential drops in the electrical 
double layer during the extreme cathodic polarization to which the 
stated condition corresponds, it would be expected that small 
inorganic cations (in this case K ) and even water dipoles could 
displace the large quaternary amine cations. 

Less conclusive experiments with tetraethyl ammonium 
bromide indicate somewhat similar behaviour but with desorption 
occurring at less cathodic potentials than with the cetyl ammonium 
bromide. This result would argue for the case of the cetyl quaternary 
amine in preference to a smaller one for dendrite inhibition. 


REFERENCES 
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CONCLUSIONS 

(1) NR^ + ions are suitable for dendritic growth inhibition 
because of the highly -ve position of the potential of zero charge. 

(2) In the absence of satisfactory results for the concen- 


tration dependence, we do not know if the adsorption effects are 
heterogeneous . 



9 


FURTHER WORK NEEDED 

The ideal dendritic inhibitor would be one in which there 
is a minimal concentration in solution for blockage of the rotation 
of spirals which is the preliminary to dendritic growth. 

We need, therefore, a study of adsorption heterogeneity to 
determine where the ions are going, and we need examination of a 
series of cations to get the law between inhibition and size. 
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WORK OF MR. R.K. SEN 

THEORETICAL CONTRIBUTIONS TO ELECTRODE PROCESSES 

GENERAL DESCRIPTION OF WORK. 

Mr. Sen was put to work on a variety of theoretical 
problems which cameup in discussions during 1970 and 1971 of 
batteries and their problems. Tn the following ( the work is 
reported . 

A summary is then given of the findings. 
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VARIATION OF THE COEFFICIENT OF FRICTION WITH 
POTENTIAL FOR A SOLID-SOLUTION CONTACT: 

A REVISED CALCULATION. 

Several authors 1 ’ 2 have pointed out that a linear relation- 
ship exists between the coefficient of friction for a wet contact and 
the associated potential difference. Bockris, Argade and Gileadi 3 
reported experiments showing a parabolic relation between friction and 
the potential difference. Bockris and Argade 4 then attempted a 
quantitative calculation of the dependence. The essential aspect of 
their model was as follows: 

The coefficient of static friction is given as: 

u f = F/R, 

where F is the frictional force, i.e., the tangential force required to 
separate the region of microcontact. Then F = OjA^, where a is the 
tensile strength of the softer of the two metals in contact and A Q is 
the area of microcontact. R is the normal reaction, given as: 

R = [w cos9 - F (A - AT] -pA + F . (A - A ), 

L a^ o ' 1 K o o el v o'’ 

where w is the weight of the slider, 0 is defined in Fig. 1, F is the 

a 

attractive force per unit area between the solids not in contact, p Q is the 
attractive surface field force per unit area in the microcontact region 
and F j is the net repulsion between the double layers per unit area. 
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However, it can be shown that if Amontons law is valid 5 , i.e., 

the coefficient of friction is independent of the geometric area of the 

contact and the weight of the slider, then 

w cos0 + F ,(A - A ) > p A + F (A - A ) . 

el o' "o o a o' 

So eq. (1) then becomes 


Vo 


q f w cos0 + F 


el 


CA 


A ) 
o' 


(3) 


To estimate F g ^ Bockris and Argade's 4 analysis involved the 
calculation of the repulsion of two plane sheets of charge, i.e., the 
other side of the charge on each double layer. This probably is an 
unacceptable simplification. In the present note, the problem is 
resolved on the following model: it is assumed that in the double 
layer concerned, specific adsorption is relatively small. In that 
case, it is possible to regard the interaction approximately as that 
between two Gouy layers*. This repulsion can be calculated by getting 
the double layer field due to one double layer at x and finding out the 
force it exerts on the diffuse layer charges of the other double layer, 
taking t0 be localized as a sheet of charge at 1 /k distance from 

the solid surface - see fig. 2a. Hence, 


el 


= ^diff C " e " KX) + K ^ e " K(X + 1/lC)) * 


M 


Under the assumptions made above Hence, 


„ . z -kx -k(x + 1/k)-, 

F el * "difrt,! 6 - e !• 


Substituting 


in eqn. (.5) yields 


q diff 


(2 H) 2 ° kT 


(4) 


( 5 ) 


* Repulsion of two Gouy layers was also calcul-ted by Frumkin 6 using a 
different approach. 
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F el = K 


en°kT 


2A 


"2 7 - 

ty 2 

kT 


exp{- kx} - exp 


- K 


1 

x + - 
Kl 


( 6 ) 


On the other hand, expressing (Op o> eq. (5), in terms of q diff yields 

4 ^ q diff| 


el 


exp{- kx} - exp 


f ll 

- K X + — 

l kJ 


(7) 


Introducing eq. (6) or (7) in eq. (3) yields eqs. (8) and (9) 


V f - 


a T A c 


W COS& + K 


en°kT 


2A 


Ze 


and 


P f = 


o .2 r -KX -k(x + 1/k),TA - A ) 

-£f [e - e o' 


¥\> 


w cose + (47Tq 2 diff /t) [e" Kx - e' K(x + 1/lc) ](A - A Q ) 


6 ~ 2 

The weight of the slider was 7.1 g, a T = 1.266 x 10 g cm ; 


A q was evaluated by Argade 6 7 and was found tfy be 7.2 x 10 


-7 17-1 

evaluated for x = 6 x 10 cm and < = j x 10 cm ; ‘Ijiff 


( 8 ) 


(9) 


Eq. (9) was 

2 

was obtained from 

the integration of the C-V plots. The results are shown in fig. 3 and 
compared with experiment. The agreement seems fair. 

Argade* s experiments further showed that when the ring was made of 
quartz a similar dependence of on potential was obtained. The explanation 
of this phenomena can also be obtained from the above theory. The quartz 
being an insulator can be approximated to a very lowly conducting semi- 
conductor having a very large Debye length (1/k^) . This can very easily 
be shown, e.g., for an intrinsic semiconductor 


n in = ( N c N V ) ^ (exp[ - E / 2kT]) 


and 



14 - 


k, = [83&i. e^/ekT]^. 

1 L m o' 1 

4 

Taking e = 4 and E^, the band gap, of 0.7 eV, k^ comes out to be 5.85 x 10 

-1 -4 

cm , i.e., 1 /k = 1.7 x 10 

Now the model is as shown in fig. 4 and following essentially 
the same principles used above, we can write 


el 


en kT 
o 


A 


Ze 


— ifi^[e~ KX - e -K( ' X * Vicpj 


kT r o‘ 


-7 -4 

Now can be evaluated for 1/K = 3 x 10 cm, 1 /k^ - 1.7 x 10 cm and 

_7 

x = 6 x 10 cm. The value of F g ^ was used to calculate from eq. 

(3). The values are plotted in curve 2 of fig. 5. The agreement again 
is fair. 


REFERENCES 

(1) J.O'M. Bockris and R. Parry-Jones, Nature, 171 , 930 (1953). 

(2) D.N. Staicopolus, J. Electrochem. Soc., 108 , 900 (1961). 

(3) J.O'M. Bockris, S.D. Argade and E. Gileadi, Electrochim. Acta, H, 

1259 (1969). 

(4) J.O'M. Bockris and S.D. Argade, J. Chem. Phys., 5<D, 1622 (1969). 

(5) J.A. Kitchener and A.P. Prosser, Proc. Roy. Soc. (London), A242 , 

403 (1957). 

(6) A. Frumkin and A. Gorodetzkaja, Acta Physiochem. USSR 9^ 327 (1938) 

(7) S.D. Argade, Thesis, University of Pennsylvania, 1968. 














® © 

© 


© 

' © © 

© 


' ~ 0* 

© © 

© 

_ © 
© ® 

© 

Solution 


Figure 4. 

■ate 








Coefficient of friction 



Fbtential(NHE)V 




15 - 


Figure 1 

Figure 2 

Figure 3 
Figure 4 

Figure 5 


ANNOTATIONS TO FIGURES 


Forces involved in friction. 

Model of interactions of the Gouy layers between two metal plates. 

Coefficient of friction versus potential for Pt on Pt in HCIO^. 

Model of two interacting Gouy layers between a metal and a 
semiconductor. 

Variation of the coefficient of friction with potential at metal- 
quartz junction. (a) Experiment; (b) theory. 
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DEVELOPMENTS OF THE BDM MODEL OF THE DOUBLE LAYER 
Abstract 

Following the publication of the Bockris-Devanathan -Muller 
model of the double layer, several discussions in support and against 
various aspects of that model have appeared. A critical discussion of 
the present status of the model is presented. 

I. INTRODUCTION 

In 1963 Bockris, Devanathan and Muller 1 (BDM) published an 

attempt modelistically to interpret several facts, at that time 

unrationalised, concerning the double layer at a mercury- solution 

_2 

interface, namely: (1) The capacity is constant at 16pF cm for the 
double layer region on Hg, whilst the cation size is varied. (2) There 

is an inflection in the experimental c l an i on -c l ine . ta ^ relation. (3) There 

is a hump on the capacity-potential curves. (4) The relation between 
coverage on the surface of Hg with organic molecules and the charge of 
the electrode is parabolic and symmetrical, if small aliphatic molecules 
are being adsorbed. 

BDM explained the independence of the double layer capacitance 
upon cation size by modifying a proposal made by Watts-Tobin, Mott and 
Parsons 2 , where it was suggested that the capacitance could be explained 
in terms of that of a low dielectric water layer; and a high dielectric 
second layer, further out towards the solution, in series with the low 
dielectric layer. The shape of the isotherm and the capacitance hump 
were explained in terms of anion repulsion among the specifically adsorbed 
molecules. A calculation of the charge at which the capacitance hump 
occurs has been shown to agree well with the available data 7 . The parabolic 
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dependence of coverage on potential for the adsorption of aliphatic 
organic molecules was interpreted in terms of the dependence of the 
strengths of attachment of the water dipoles to the electrode, the dipoles 
having a two state configuration. It is the turning round of the water 
molecules which causes the maximum in adsorption of an organic. The 
presence of quadratic terms in the field, arising from the polarisability 
term, was shown to be numerically insufficient to explain the observed 
parabolic 0-V shape. Bockris, Gileadi and Muller 3 (BDM) published a 
better approximation of the BDM model for organic adsorption, in which 
both solvent solute dipoles interact with the field. Coverage-charge 
relations of numerous shapes were deduced. However, the lateral inter- 
actions between adsorbed organic particles and between organic and water 
molecules, were still neglected. 

Thus, the BDM model features potential-dependent water 

molecules as a primary constituent of the double layer; differs from a 

previous model 2 in holding the distance of closest approach of non- 

specifically-adsorbed ions as outside this water layer; and attributes 
the capacitance hump to anion repulsion rather than to a contribution 

to the capacitance due to water molecules 6 .* 

Support for the BDM model in terms of the isotherm 7 was also given 
by Wrob Iowa and Muller 8 , who showed that the model's neglect of imaging in 

* A water capacitance hump is also implicit in the BDM model. Whether one 
observes a capacitance hump due to anion repulsion, or to water orientations 6 , 
turns out 1 to depend on the water-water interaction energy, which has been 
neglected in the approximations of Mott and Watts -Tobin 6 and of MacDonald 
and Barlow 23 . In BDM's evaluation, these interactions are shown to make 
water give a capacitance too large to be observed. Anion humps are 
qualitatively consistent with the qualitative facts of hump-dependence 
of anions. Whether capacitance humps are always anionic will depend on 
solvent as well as anion. 
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the solution gave a better agreement with experiment than the consequences 
of the assumption that there is perfect electrostatic imaging at both the 
metal surface and the OHP (which involves the assumption that there is a 
sharp electrostatic discontinuity between the double layer and solution 
regions). Gileadi 9 has shown that the terms (of water-electrode interaction) 
used in the BDM model for organic adsorption, have a much greater numerical 
weight (in the potential range accessible experimentally) , than terms 
representing the change of energy of the double layer which arises in the 
classical thermodynamic approach 18 to the coverage- charge relation of 
organic adsorption (cf. the earlier discussion by Wroblowa and Green 17 ) 

Gileadi and Stoner 10 have developed the BDM model, with respect to the 
dipole contributions to the double layer potential, as an interpretation of ano- 
malous current potential relations. Reddy 11 has reaffirmed the model, with 
special reference to the water contributions. Barradas 12 has argued that 
BGM 3 model represents consistently some phenomenology of organic adsorption. 
Bonciocat 13 has shown that the BDM isotherm gives rise not only to the 
hump, but also to the minimum in capacitance observed at potentials positive 
to the hump, and that the isotherm rationalizes the temperature dependence 

Wan 

of the hump. Isotherms which take into account multiple imaging, such as i 

that of Levine, Bell and Calvert, or MacDonald and Barlow, do not do this. 

Levine, Mingus and Bell 14 , and MacDonald and Barlow 15 , had criticised 
' • ■ Ij ' 

the BDM model because of its lack of multiple imaging (cf . Wroblowa and 

Muller 8 ). An alternative theory 6 , in terms of capacitance humps in 

terms of water rotation capacitance, competes with the BDM model. V 

Damaskin 16 has suggested that consequences of the BDM model of organic 

adsorption are inconsistent with thermodynamic reasoning. He has suggested | 

ft 
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that a neglect of the energy of charging the double layer is important 
(but cf. Gileadi 9 , Wroblowa and Green 17 ). The neglect of electrode- 
organic interactions in the original BDM has been criticised by Frumkin 18 , 
but this objection is no longer relevant because a second approximation of 
BDM (that of GBM) took this into account. The neglect of a polarizability 
term in %aF has been criticised also by Damaskin 16 '. Mohilner 19 has 
suggested that the dipole potential arising from the BDM model is too 
large. 


11 . DEGREE OF CONSISTENCE OF THE BDM MODEL 

1. The Ionic Adsorption Isotherm 

The BDM isotherm can be represented in the form: 

= -AG°/RT + Zn— + Aq., + B0 3/2 . (1) 

1-0 a n M 

o 

This isotherm is obtained using single imaging for the 
evaluation of the lateral interaction terms. Levine, Bell and Mingus 14 
and MacDonald and Barlow 15 preferred multiple imaging. The multiple 
imaging results in a higher screening effect, decreasing the lateral 

, 3/2 

interaction energy, which then becomes proportional to 8 instead of 0 
obtained by BDM. Thus, the Levine, Bell and Calvert 28 isotherm is of the 
form: 

to* S' = “ ^§f + + Ac *m " C0, (2) 

(l-p6) p K1 a o ™ 

where p = r . /r . (assumed approx. 1 in BDM). 
r anion water rr 

Which of the two approximations for imaging agrees more closely 

with the physical situation? The tests that can be carried out are (1) 

plotting the log a + (l-0)/0q M vs 0 3 ^ 2 and log a. + (l-p0) P /0q M vs 0, to examine 

the linearity as well as the slope from experimental data, and (2) by 

8 

examining the Essin-Markhov effect. Wroblowa and Muller made these tests 
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and their results can be summarised as (1) the slopes B and C of 

equation (1) and (2) when compared with experiment from the above 

mentioned plots, show that the predictions of B from BDM model agree 

within about 15-20% of the experimental results, whereas the C values 

predicted by the LBC isotherm, taking p = 1, exceed the experimental 

ones by four to twelve times for the anions Cl , CN~, C10~ 4 ,l, 

CIO^ , BrO^ and CNS , and (2) for p =i the constancy of 

dqCA/d q M at covera S es usually encountered in ionic adsorption is 
a ± 

not obtained from the LBC isotherm, whereas the same quantity from the 
BDM model does almost stay constant. 

The single-imaging approximation gains in these tests. 

Single and infinite imaging are limiting cases of a more complicated 
situation which neither equation (1) nor equation (2) represents 
accurately. It is more difficult to think, however, that conductive 
imaging (the assumption of LBC) in the relatively poor conductor, an 
ionic solution, will occur to a significant extent. MacDonald and 
Barlow 15 pointed out that, even in the metal, imaging at molecular 
distances is not perfect. Due to the lack of sharpness of the dielectric 
boundary between the OHP and the rest of the solution, dielectric imaging 
beyond the OHP will be far from perfect, and that is probably why single- 
imaging gives better agreement with experiment. An estimate of 
partial imaging, on the solution side, is desirable. Some attempt has 
been made in this direction by Buff and Goel 29 . In concurrence with the 
present discussions, their results show that the more gradual is the 
change in the dielectric constant the less pronounced is the degree of 
imaging in the solution. 
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Recently, Levine 30 has recalculated the results obtained 

by Wroblowa and Muller 8 from his multiple imaging model, but using 

a value of p = 2, for his isotherm. The term p represents the radius 

ratio of an ion to water and it is nearer to 1 than 2 for most anions 

concerned. The use of p as an adjustable parameter and its substitution 

in LBC decreases the discrepancies with experiment which this model has 

for p = 1, so far as the slope and the Essin-Markhov effect goes. 

(It is clear that all isotherms should be formulated with a p value of 
2 2 

r l ^ r H (P ' ^ e » t ^ iere ^ ore » compare the single and multiple imaging 

models in respect to the capacitance hump phenomena. In respect to 
the BDM model, the vs curve (and thus the hump, with anion 
dependence) is predicted 7 , and recently it has been shown 13 that both the 
capacitance hump and the subsequent minimum are predicted by the BDM 
theory. Let us see whether the prediction of humps from an LBC isotherm 
is also consistent with experiment. The LBC isotherm can be written 
(p = 2) as: 


Hence : 


A’q M = ln0 - ln(l-2&) - lna + - C& + constant. 


dc k 


a i n M _ 1 4(1-29) 

d0 ' 0 ( 1 - 2 S ) 2 


(3) 


So: 


1 + 


(0 ( 1 - 20 ) 


- C 


(4) 


d9 


A 


I + 4 

0 1-20 


- C 


( 5 ) 


Let x Q + x _ 2 q - c 
Then: 


d0 
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M 


Al 

x 
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Now: 


d 2 9 
d< >M 2 


2 


d 

de- 


de 


dq : 


M 


(A*) 2 dx 

3 ' d8 ' 

x 


Let us now evaluate 4^- . 

do 


So: 


dx 

de 


d 

de 


F* 


8 

d-26) 


(1-26) ‘ C J 

46 2 t 46 -1 

2 2 
6 (1-26) 


( 6 ) 


Thus: 


d 2 6 _ (A') 2 

d V 3 


46 +46-1 


l 0 2 d - 26) 2 J 


d 2 e 

For the inflection point dq M 2 


= 0. 


(7) 


Thus: 


(A') 2 

2 ' 
46 -46-1 

3 

X 

. e 2 (i - 28 ) 2 


( 8 ) 


A' t 0. Thus: 46 2 + 46 - 1 = 0 . 

Or: e x = |(/2 - 1) or 6 2 = - |(/2 + 1) . (9) 

A negative value of 6 has no physical significance. Thus, the inflection 
point is given by the root 6^. Hence (cf. relation (9)), the inflection 
in LBC is independent of any property of the ion, or the value of q^. Thus, 
for ions to which the LBC isotherm is intended, the plot of q CA vs q^ the 
inflection should occur at the same point for every anion. This is not 
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the case. The plot of 1i n f lection vs \ ump woul <i not give a straight line 
with a slope of unity as is obtained from the BDM theory. 

Thus, the multiple imaging approach of Levine et al 29 does give 
less consistence with experiment (qualitatively as well as quantitatively) 
in its application to predict the inflection on the plot of q^ A vs q^. 


2. The capacitance hump 

The most discussed aspect of the recent double layer theory has 
been the rationalization of the hump observed in capacity-potential relation. 
The BDM model formulates in terms of single imaging an idea suggested 
qualitatively by Devanathan 20 and supported by Watanabe 21 . In BDM, it is 
the contribution of the repulsive forces to the standard free energy of 


adsorption which provides an inflection on the q 


CA 


q^ curve, and thus gives 


a hump in the C - q M curves. The BDM ionic adsorption isotherm for low 

3/2 


coverage can be written as: 

‘ 0 


In 


1-0 


= constant + In a. + Aq 


M 


B0' 


( 10 ) 


where constant = 


AG° 

RT " ma o 


The condition for the inflection point in the 0 = f(q M ) curve can be 
written as: 


2 - 40 - |B0 3/2 + B0 5/2 = 0, 


(ID 


where: 


B = 


7/2 2 2 

. 7T ' ,e r. 
o 1 


4ekT 


Snax^A 


3/2 


where the terms have their usual meaning. 
This equation has two roots 13 : 
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0j_ = (4/3B) 2/3 and 0 2 = 1/2 . 

corresponds to a local maximum, i.e., the hump in the C - q M curve, 
and 0 2 is the local minimum. 

On the other hand, considering the root of the LBC isotherm, 
it can be shown 13 that it does represent the maximum. But, since the 
root 0^ of the LBC isotherm does not depend on the nature of the anion, the 
capacitance hump in this approach is independent of the nature of the anions 
which is at variance with experimental observations 21 . The LBC isotherm 
furthermore does not predict any minimum after the hump, whereas the BDM 
model does. 

Others 22123 consider the hump is due to a solvent effect. 

There are qualitative objections to the solvent reorientation model for 
the hump. (i) For common anions except F , the region of charge near to 
the hump sees the onset of anion adsorption 24 ; (ii) the magnitude and 
location of the hump varies with the anion 25 ; (iii) in some polar non^- 
aqueous solvents, no hump is observed 25 ; (iv) if solvent reorientation did 
give rise to a maximum in local permittivity, and hence the hump, it would 
do so for aqueous solutions on the cathodic side of the pzc 24 . 

In favour of a solvent reorientation theory of the hump is that 

humps occur for F . There is considerable evidence 26 that the F is not 

specifically adsorbed. However, this conclusion relies on the applicability 

of the Gouy-Chapman diffuse layer theory, the assumption of zero specific 

adsorption of water, etc. It is possible that F merely differs in 

degree of adsorption from other halide ions. The evidence is strong that 

specifically adsorbed anions play a role in causing the hump. The BDM model 
quantitatively predicts the phenomenology. 
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Payne 27 found in N-methyl formamide system two humps, one 
on the cathodic side and the other on the anodic side of the e.c.m. Both 
BDM and solvent reorientation theories find difficulty in explaining both 
humps. However, if solvent-solvent repulsion is relatively small, BDM 
also predicts a solvent capacitance hump as does the interaction-free, 
anticedent model of Mott and Watts-Tobin 6 . In two hump systems, one 
hump may be due to anion-repulsion, the other to dipole orientation. 

3. Organic Adsorption 

The BDM competition theory of organic adsorption is a development 
of the water competition model proposed by Butler 4 . It will not be 
applicable to potential regions in which there is significant chemisorbed 
hydrogen or oxygen on the surface, nor to situations in which there is 
dissociation of the organic molecule on adsorption. The original version 
of the model 1 will only be applicable to systems in which there is no 
electrical interaction of the double layer with dipole moments of the 
organic molecules. This assumption applies to a situation in which the 
organic has its polar groups in the diffuse layer, where the field is 
negligible compared with that on the water dipole, which is inside the 
compact layer. In the BGM 3 approximation, various orientations of the 
adsorbed species were considered and then interactions with the electric 
field explicitly calculated. The model was able to interpet most of the 
shapes of the 0-q M curves obtained. 

(i) Is the BDM model inconsistent with thermodynamic reasoning ? 

Let us take a general isotherm, congruent with respect to charge 


of the form: 
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G(q)c = f(D (12) 

where £(r) is a function of the surface excess, T. 

Thermodynamic treatment 16 yields an equation of the form: 


G q exp -A 


-l CK-c’) 


— q + Axdq + qE° rg , 


where G q is the value of G at q = 0, A = RT F, K q is the capacity of the 
double layer when the coverage of the organic, i.e., 0=0 and c' is the 
capacity when 0 = 1, q is the charge. Ay is the surface potential and 
E° rg is the potential drop due to the oriented adsorption of the organic 
dipoles . 

Similarly, if we take the BDM organic adsorption isotherm which 
is: 

C -AG° /RT nAG° /RT - (Ec/kT) , , . . REC) .... 

9 org e org g c' e 1 } !(4TTyq/e) - , (14) 

1-0 C W 

we can rewrite it in the form: 


' \ 

G(q) = G exp -nR(yX - REC)/kT (15) 

0 l 

47rq M 

where X = electric field strength = — — , 
and R is given as: 

R = tanh (yX - REC)/kT 

where G(q) and G q have the meaning of equation (13), n is the number of 
water molecules displaced during the adsorption of a single organic 
molecule, y is the dipole moment of water, X is the electric field acting 
in water, E is the lateral interaction energy between water molecules, and 
C is the effective coordination number of the water molecules at the electrode. 
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Equations (13) and (15) differ importantly. 


The terms 


.-1 (K o - c ’> .2 

1 2K x ' q 

o 


and 


qE° r ®j are absent in (15) . Since 


V. , * / 

equation (13) is derived 6 from a charge-congruent isotherm, it should 
be thermodynamically consistent solution of the problem. As the BDM 
isotherm, represented by equation (15), omits two terms of the isotherm 
of equation (13), it has been argued 16 that the BDM isotherm is incon- 
sistent with thermodynamics. 

In the derivation of equation (13), given by Damaskin 16 , 
one involves the equation: 


E = E 0 (l - 0) + E0 (16) 

where E is the total potential*, E' the potential at 0 = 1, and 
E q is the potential when 0=0. This equation is based on the 
modelistic assumption that, at a given charge q, there is a linear 
variation between E and fe. This will only be true if the free energy of 
adsorption is not coverage dependent, which is not generally true. 

This fact has been pointed out by Delahay 31 . Thus, equation (16) is not 
purely thermodynamic. It represents a model. Consequently, equation 
(13) cannot be purely thermodynamic. Thus, both equations (13) and (15) 
involve modelistic assumptions and, as these differ, their consequences 
would not be expected to he the same. 

* E in the Russian literature is defined as the total potential. 

It is implicitly assumed that the rest of the potential in a cell 
remains constant so that change in E reflects the change in the potential 
at an electrode solution interface. 
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An argument of a similar nature arises because equation (13) 
is obtained from equation (12), which is an adsorption isotherm based on 
a model (essentially, electrode-water and electrode-organic interactions 
are considered), and equation (15) arises from a different model (essentially, 
interaction of water molecules with the electrode and each other) . The 
first and third term of equation (12) correspond to a change in the double 
layer energy due to organic adsorption, and the reorientation energy of 
the organic dipole, respectively. However, the BDM model involves an 
initial assumption of a condition such that these energies are smaller 
than the electrostatic interaction of water (and, in BGM, the organic), 
with the electrode. It may be right to argue that the neglect of change 
of energy in the double layer due to change of capacitance on organic 
adsorption is a poor approximation (see next section) . However, it is 
not rational to claim, as did Damaskin 16 , an inconsistence with thermo- 
dynamic principles because another model (which has been formulated in a 
way which involved thermodynamic equations) predicts some other behaviour. 

(ii) Neglect of the change of the energy of the double layer upon 
adsorption of the organic : 

It has been argued 12 that this approximation in the BDM theory 
could lead to important discrepancies with experiment. Recently, Gileadi 9 
evaluated the relative contribution of the term corresponding to the energy 
change due to water dipole reorientation, and the term due to the energy 
change of the ionic double layer due to organic adsorption. The results 
are summarised in T able I. 
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TABLE I 


Relative contributions of the Frumkiri term and the Bockris 


Potential in volts 

et al 
0.5 

term 

1 

2 

5 

10 

Term corresponding to 
energy change in ionic 
double layer .% 

7.6 

14 

25 

45 

62 

Water dipole re-orienta- 
tion term % 

92.4 

86 

75 

55 

38 


Table I shows the relative contributions of the two terms 
to the total variation of the adsorption equilibrium consistent with 
potential at different values of E. Although the term corresponding 
to the change of energy of the double layer is not negligible (as 
assumed in BDM), its contribution will be relatively small over the 
potential range accessible experiment (M volt). Thus, this objection 16 
is not a weighty one. 

(iii) Is the Electrode Charge predicted for the Adsorption Maximum 
BDM (or BGM) inconsistent with experiment? 

It has been inferred 16 that, since in the BDM theory the 
maximum of organic adsorption occurs at the point where the 

maximum of organic adsorption should occur only at the pzc. Since 

experimentally this is not observed, the argument was that the BDM model 
does not predict the experimental results well. This criticism is due 
to a misreading of the original papers. In BDM, although it is argued 
that the maximum of organic adsorption occurs at the point where = N^, 
it is pointed out that this situation should not arise at the pzc, because 
of the stronger binding of the water molecules in one orientation than 
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in the other.* Bockria and Swinkels 32 , using the BDM model, showed 
that the maximum of organic adsorption will occur when: 


(AG°+) - (AG°I) 
L C 


8Try 

e 


Snax 


(16) 


where AG^i and AG° + are the free energies of adsorption of water in up 
and down position. 

This is expected for a model in which there is no electrical 

interaction of the organic molecule with the electrode. For the more 

general (BGM) model 3 , in which such an interaction is accounted for, the 
adsorption maximum will depend not only on the water-electrode, but also 

organic-electrode interactions. It will, finally, depend on all 

interactions (e.g., also organic-organic and water-organic), although 

others may not have a great effect (see Section V) on the adsorption of 

aliphatic molecules. 


(iv) Is there any relation between the hump on the capacitance vs electrode 

charge curve and the maximum of organic adsorption? 

It has been argued 16 that there is a relation predicted by 

the BDM model between the potential of the hump of the capacitance curve 

and the potential of maximum organic adsorption. The BDM model predicts 

that the peak of organic adsorption occurs at the point where = N^. 

Since the capacitance hump (if it is attributed to the solvent) would 

also occur at the point where = N^, the two potentials ought to be 

the same. Experimental results do not show such a correlation. Hence, 

the BDM model is inconsistent with experiment. 

* The difference may arise, inter alia, from the lack of symmetry and hence 
water-metal image energy, in the two positions. 
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However, this, argument misrepresents the BDM model: the 
capacitance hump in this model is not generally due to solvent reorientation, 
but anion repulsion (Section II, 1). There is absolutely no relation, 
in the BDM model, between the potential or the electrode charge, at which 
there is an adsorption maximum, and that at which a capacitance maximum 
occurs - in agreement with experimental observations. 


(v) Is the BDM model neglect of organic-organic and organic-water 
interactions significant? 

The neglect of the organic-electrode interactions (one of many 
approximations in the BDM model) was revised in BGM. Bockris and Swinkels 32 
made an estimate of the organic-organic interactions, but they have not 
yet been introduced into the general equations of improvements in the 
treatment of the BDM model, largely owing to algebraic difficulties. 

The BDM model has been criticised for neglect of organic-organic 
interactions. Let us, therefore, make a crude estimate of the importance 
of organic-organic interactions in organic adsorption. We will assume for 
simplicity that there are only two positions of the organic molecule. The 
interaction of a pair of these, when they are pointing in the same direction, 
is : 


- V 

“12 3 

er 


(17) 


Due to the effect of the value of e, E ^ for water is about 7 times larger 
than that of the butanol-butanol interaction. Even for 0 = 0.5, where 
the distance between water molecules and the organic molecules are 
approximately the same, i.e., r is the same, the organic-organic interaction 
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energies are one-tenth of the water -water interactions. Only when the 
coverage by the organic approaches 0.8 does the organic-organic inter- 
actions become of the comparable order of magnitude as that of the water- 
water interactions. Thus, neglect of the organic-organic interactions 
in studying organic adsorption seems a reasonable compromise with the 
complexity of the expressions relating coverage to electrode charge. 

However, organic-organic interactions at higher coverage may dominate 
the interaction situation where the organic dipole lies on the electrode 
(see Fig. 1). In such a case, the objection is valid. 

In a recent article Dutkewicz, Garnish and Parsons 34 showed that 

for 1:4 butane diol and 1:4 but-2-yne diol the charge dependence of the 

free-energy of adsorption does not give the same result as obtained from 

the Bockris, Devanathan and Muller theory. They have interpreted the 

discrepancy as arising due to the neglect of the water-organic interactions. 

It is clear that organic-water interactions should be accounted for in 

principle, and that they will become significant if the dipole of the 

organic sits in a position where it is exposed to a field of about the 

same magnitude as that of the water. The BDM formulation would apply when 

the organic dipole is in the diffuse layer and the approximation gets 

better in the organic dipole away from the water. The discrepancy would 
be particularly great if the organic molecules are lying flat on the 

electrode surface. 

However, physically when the organic polar groups are in the 
high dielectric constant region, whereas the water dipoles are in the low 
dielectric constant region, there seems to be little reason why the 
interaction should become very important. 
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(vi) Neglect of the polarizability term in electrode molecule interaction : 

_2 

It is trivial to show that, for (q M )<10y coul cm : 

j(Aa)x 2 «(Ay)x 

4nq M 

where x = — — and e is taken as 6. This conclusion assumes the two-state 
water dipole model, and in both, the dipole is perpendicular to the 
electrode. 

-2 

At q M >10y coul cm , the polarizability term will become 

increasingly significant and should be included in a further appoximation 

for the model. However, it is important to note that the small value of 
2 12 

l/2ay (or ^CV ) compared with yx> at X ^ means that change in direction 
of the converage- charge relation cannot arise (as in the former models) , 
by virtue of the quadratic term . It is therefore necessary to have some 
change in position of an electrically interacting entity in the double 
layer, near the adsorption maximum. This was the origin of BDM's suggestion 
that this change is in the orientation of the water dipole.* 


4. Surface Potential : 

Mohilner 19 stated that, assuming full coverage of the electrode 


* The two state model assumed for water is, of course, also an approximation 
(cf. MacDonald and Barlow, who discussed a multi-state arrangement). 

The effect of several molecular orientations in respect to sample organic 
molecules is formulated in BGM. The symmetry of the 0-q^ relation often 

observed does, support the two-state model for water in BDM. Howev er, other 
possibilities (which would reduce the dominance of the water effects) are: 

(a) The introduction of a degree of dimerization; (b) A two state model in 
which the effective dipole moment is reduced due to orientation other than 
vertical. 
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with water molecules, the value of the surface potential Ay from the 
BDM model is unrealistically high. 

In the BDM theory (and for an organic of zero dipole moment) : 


A X = 


. 4iryN n 


tanh 


UX REc 
kT “ kT 


where R = tanh 


UX nl£ 
kT kT 


(18) 


Moreover, to eliminate the water hump in the capacitance 
dy 

potential curve, A dipole „ . -2 

? » _c 0.01 when q M = 20y coul cm . 


i.e. , 


; M 


8 ><dlp 1 

q M C dip ~ 


(1-R)' 


< 0.01 . 


s kT _(Ecl (1-R ) 

kT 


(19) 


In judging what parameters should be used to test this equation, 
there are two constraints. The change of Ax with potential must not be too 
high, because if it is, it will not be consistent with the observed linearity 
of the Tafel line. 

N0W ‘ 4ttu h q N t 

AX = 7 T tanh <0.1 volt 


UX REc 
kT " kT 


and 


I 

q M 


C dip 


8TT 2 yN T 

e 2 kT 


1-R 


1 + ff-Cl-R 2 ) 


« 0.01 . 


Evaluating these two equations with the above constraint yields 
14 2 E 

Ay - 100 mvftr N^, - 10 molecules/cm for pjr = 0.5 (the value assumed by 

BDM) . The problem is that the value of N,p thus required is significantly 

less than the close packing value for the water molecules. To obtain a 

total coverage of only c.10 1 ^ molecules/cm" 2 , r u n * 3.9% instead of 1.4$. 

**2 • 
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How may we reduce this difficulty? Examination of the above 

E 

equation shows that, as we increase the value of pjr , the value of 


C,. 

dip 


1 E 

, and also the value of Ay, also decreases. For -1, the 


value of Ay - 100 mv is obtained for N T = 5 x 10 14 molecules/cm 2 , 
which is the close packing value if we assume the radius of water 
to be - 2 X. The problem, therefore, is physically to rationalize 

E 

the value of pjr - 1. Let us assume that the lateral interaction 
energy is given by dipole-dipole forces plus the quadrupole-quadrupole 
forces . Thus : 

2 2 

E = E ... + E„ = H- + ? 


yy 


QQ 


R" 


R" 


Evaluating this for interaction between two water molecules separated 
by 2$, yields E - 0.5 Kcal/mole and Eqq = 0.08 Kcal/mole. So, E = 

E 

0.58 Kcal/mole. Hence, pjr = 0.966 which is near to the value needed 
to give reasonable agreement for the close packing value of N^,. 

Thus, a reasonable value for Ay does arise on a BDM model if a better 
approximation is used for E, the water-water interaction energy. 

III. PRESENT POSITION OP THE BDM MODEL 
The BDM model gives a prediction of humps with those observed: 
its prediction of linearity in log a + - log — - versus 0 2 is consistent 
with experiment. Some degree of multiple imaging should be introduced, 
but it must be less than that of the multiple imaging models to retain the 
consistencies of the single imaging model. The BDM model is consistent 
with a solvent hump, if the lateral interaction among the solvent molecules 
is low enough. 
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For organic adsorption, some of the criticisms published 
involved misunderstanding of the model, e.g., that about the position 
of the hump, the peak of organic adsorption, and that which neglected 
the BGM organic interactions. The most serious criticism, that the 
model was in disagreement with thermodynamics, is not correct: BDM is 
not consistent with another model, which has used thermodynamics in its 
formulation. The neglect of the ionic double layer energy seems a 
reasonable approximation over the potential range accessible. 

The three principal contributions of the BDM model are 
the consequences of interactions among water in a two state model; 
the demonstration that such a model rationalizes the 0-q^ parabolas 
for aliphatic adsorption (and that models which rely on quadratic terms 
in V or q will not do this); the position of the outer Helmholtz plane 
as being outside this layer; and the demonstration of a reasonable 
degree of consistence with experiment for single imaging. These 
contributions seem not to be challenged by valid criticism. The most 
important directions in which the model should be developed is in respect 
to a better evaluation of the water-water interactions, with an examina- 
tion of the effects of various other configurations for water, the 
inclusion of organic-water and organic-organic interactions, and perhaps 
an attempt at a realistic multiple imaging calculation. The last 
modification seems the least needed (in respect to consistency of ionic 
isotherms with experiment) and the most difficult 33 . 

The continued neglect of contributions to the double layer 
properties of the potential dependent water dipoles, except by MacDonald 
and Barlow 23 , seems the major anomaly of the last decade of double layer 


research. 
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IV. GENEALOGY OF THE BDM MODEL 

Lange and Miscenko 41 suggested the presence of a surface- 
potential, Ax, in double layer studies. A formulation of the water 
layer at an electrode solution interface in terms of a dipole layer, 
was first formulated quantitatively by Bockris and Potter 35 in connec- 
tion with the interpretations of pH effects in the hydrogen evolution 
on Ni. They were the first to introduce a field dependent contribution 
to an electrode potential from the oriented water dipoles. The water 
dipole model was used by Bockris and Conway 36 (1958) to examine double 
layer properties with respect to the dependence on frequency of the 
double layer capacity. Mott and Watts-Tobin 6 (1961) picked up the 
concept of Bockris and Conway 6 (1958) ("waggling" dipoles), and suggested 

a two-state model for the water-electrode interaction to predict the 
hump on the capacity potential relations. Correspondingly, MacDonald 37 

(1959) suggested a treatment of the water layer, assuming all orientations 

of the water molecules to be equally probable, to explain the hump on the 

capacity-potential relationship. MacDonald 38 (1960) and MacDonald and 

Barlow 23 (1962) qualified the model of the water layer proposed by 

MacDonald 37 (1959). However, both Mott and Watts-Tobin 23 (1962) omitted 

the water-water repulsion, i.e., treated the layer as if it were a gas 

layer. 

BDM picked up the two state model proposed by Mott and Watts- 
Tobin 6 (1961) and introduced dipole-dipole repulsion energy. They 
showed that a two-state model was needed to give the parabola observed 
in 0 q “ Taking the repulsion energy into account helped BDM to 

observe the fact that dipole reorientation may not make significant 
contributions to the hump, because the dipole capacity was too large to 
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have a significant effect. However, BDM showed the importance of 
the dipole reorientation (reduced in magnitude by water-water inter- 
actions) in calculating the surface potential as well as in organic 
adsorption . 

BDM showed that the hydration shells of the cations must 
be outside the water layer. A layer of oriented water dipoles lies 
between the electrode and the hydrated cations. (Cf. Watts-Tobin, 

Mott and Parsons 2 ). 

BDM's proposal that ionic adsorption isotherms were better 
calculated on a single imaging theory was independent of Levine, Bell 
and Calvert 28 (1962) who had developed the multiple imaging isotherm. 
Another approach, more or less similar to the multiple imaging approach 
of Levine, Bell and Calvert 28 , but using a hexagonal lattice model was 
made by Barlow and MacDonald 42 (1964), which yielded results very much 
similar to the results obtained by Levine, Bell and Calvert 28 . 

The approach to explain the capacitance hump in the BDM model 
had an antecedent in the work of Devanathan 40 (1954), who had shown that 


1 

C 


“m- 


^M-2 Sl-l 


dq, 


CA 


dq. 


M 


However, Devanathan did not relate q CA to as was done, and compared 
with experiment, by BDM. 
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ANNOTATION FOR FIGURE 

Figure 1. Adsorption of organics on electrode: 

(1) Up position - organic dipole in diffuse layer. 

(2) Flat position - organic dipole on electrode surface. 
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QUANTUM MECHANICAL FORMULATION OF ELECTRON TRANSFER RATES 


Formulation of electron transfer rates in solution and at 
electrodes has for long been carried out in terms of continuum theory 
of dielectrics 1 ’ 2 . Some displacement, only recently defined, in the 
solvent causes the potential energy of the ion concerned to change to 
a value suitable for the acceptance of an electron, in a radiationless 
transfer process. The displacement of the solvent has been taken as 
simple harmonic, and a detailed model has been suggested 3 . A classical 
libratory movement of the solvent interacts electrostatically with the 
ion. Certain long range fluctuations in this energy, rather than 
the thermally activated vibration-rotation levels involving short range 
interaction in the ion solvent complex, give rise to the energy states 
in the ions to and from which electrons transfer. 

Two reasons have been advanced 3 for the concept that 
thermally activated vibration-rotation levels usually considered in 
solution reaction do not take part in an electrode reaction. 

First, in electrochemical reactions, variation of electrode 
potential is brought about by effecting a change of excess charge on the 
electrode from an outside source, thus causing the Fermi level to change 
by up to leV. Regarding the H^0 + levels as only those of vibration, 
there would be a gap of about 0.5 eV between levels available for the 
acceptance of electrons. The rate-potential relation is, however, a 
completely smooth one. Hence, available levels for the relevant quantum 
transition must be produced by some other means. 2 This argument, 
however, neglects rotational levels between the vibrational levels in 

* 4 * 

HgO (ref. 4), and the considerable line broadening in the spectrum of 
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HgO + caused by the introduction of solvent 5 . Schiffer and Homig 6 
interpreted line broadening of the 0-H bond, and enhanced Fermi resonance 
in water, by assuming a classical distribution for free water molecules 
giving rise to a Boltzmann distribution for the energy of the 0-H bond. 
Similar arguments may be applied to the 0-H + bond in water as well. 

Second, it has been suggested 7 that there is an insufficient 
number of thermal states to sustain charge transfer reactions . 

The following analysis shows, however, that this supposition 
is not supported by considerations of the Boltzmann distribution in 
solution to experimental activation energies. Thus, the rate of the 
proton discharge reaction can be written as 


i = F. 
= F. 


S.C ..a iS + /R .a- AH +/ RT 
if. e AS */ R N* . 


where is the number of H 0 + ions having sufficiently energetic 

H 3° 

rotation-vibration levels for radiationless transfer. For the proton 

-11 -2 

discharge reaction on Hg from acid solution, i - 10 amp cm , C ~ 

13 -2 ° H 3° 

10 1 molecules cm , and AS^ - - 10 e.u. (estimated by Temkin 8 

-3 -2 

and by Bockris and Sen, unpublished). Thus, N - 10 molecule cm . 

-i 

For AH* = 20 kcal mol (Conway 9 ), the number of particles in the second 

-2 -2 

vibrationally excited levels is about 2 x 10 molecules cm - sufficient 
to sustain the reaction. 

Thus, participation of thermally activated states in electrode 
reactions is likely 10-12 . The fact that solvent, fluctuation is unlikely 
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to be a dominating cause of activation in electron transfer reactions 
is shown by the following. (1) If one compares the "experimental" 
standard free energy of activation for a large number of redox reactions 
in solution with what is calculated on the assumption that solvent 
fluctuation is the preferred mechanism, there is little correspondence 
between the two figures 13 . (2) In electrode reactions, the logarithm 
of the rate is proportional to the overpotential applied. This relation 
does not follow from the theory of solvent fluctuations, which indicates 
a continuous curved relation between log i and overpotential 13 . 

(3) The polaron theory has been used to estimate the electrostatic 
fluctuation energy. Polaron theory is applicable in solid lattices, 
and its applicability to liquids is questionable 12 ’ 14 , and involves 
correction factors, knowledge of which is not readily available 2 . 

(4) Direct evidence for the partial nature of the contribution of the 
continuum energy considerations to charges in solution has been offered 

by Fueki et al 7 . For energy changes of hydrated electrons, approximately 

half comes from the contimuum electrostatic energy, the rest from the 
first layer of water molecules around the cavity of the electron, a 

contribution entirely neglected in the electrostatic fluctuation model. 

(5) The reorganization energy in the solvent fluctuation models has 
been given by the Born-Landau equation. The derivation of this equation 
involves the electronic adiabatic approximation, but Jortner 15 points 
out that such an approximation on polar liquids, where the binding energy 
of the solvated electron (about 1.7 eV) and that of the medium electron 
(about 4eV) are of the same order of magnitude, is not applicable. 
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The above discrepancies do not mean that some contribution 
from electrostatic solvent fluctuation does not exist. It is necessary 
to compare the probability of having electrostatic fluctuations of suitable 
energy with that of activated vibrational levels of the same energy. 

Thus, we consider an ion surrounded by N dipoles. The probability of 
the dipoles having an angle between 0 and d0 with the ion at a distance 
of R and R + dR is given as 


2r, 


P. = 
1 


"T exp 

'cd 


ey„ ff co S q 

R aR si«0d0 


EgR^kT 


TT 

r 

' °° 3 

ep eff COs6 


3 P 

2 


2r 

e„R kT 

or. (D i 

i 

l s J 


R dR sin0d0 


where r is the radius of a water molecule, r. that of an ion plus the 

diameter of a water molecule, and y g ££> following Kirkwood, is the 

effective dipole moment of the water dipoles caused by the presence 

of interactions among the water molecules themselves. Using classical 

statistical mechanics, the mean square deviation of the energy due to 
2 

fluctuation, a , can be evaluated to give 


22 2 -24 

0 = (E Z ) - (Ej = 0.01 X 10 . 

2 2 

Since a is small compared with (E) , the probability of having a 

fluctuation of 20 kcal mol ^ above the average energy can be calculated 

.-41 


using a Gaussian distribution. It is 10 


The corresponding 


probability of having an activated vibrational state of 20 kcal mol 1 is 

Mf 15 . 
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Thus the widely acceptance of the solvent fluctuation model 
for charge transfer reactions should be modified. Thermally activated 
states, as well as electrostatic fluctuational contributions should be 
counted. For the slower reactions, the solvent fluctuations do not seem 
to contribute significantly to the activation process. 
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MECHANISM OF ACTIVATION OF ELECTRON TRANSFER IN SOLUTION 
The velocity of electron transfer between a redox couple 1 2 
or between an electronic conductor and an ion in solution 4 has been 
discussed in terms of the Landau-Born 5 expression for the "electron 
solva^on energy" in solution (the continuum solvent activation model). 

It is supposed 6 that the charged particles receive energy from the 
classical librations of water molecules in a region of a few hundred 
X around the ion, in both the initial and final states. The energy 
transferred from these librators creates a situation in which a radiation- 
less electron transfer can occur. 

Acceptance of this model leads to several discrepancies with 

experiment. (1) Experimentally, the relation between logarithm of the 

rate of a redox reaction and its overpotential is linear over several 

decades 7 ; this fact is not predicted, even for small regions of potential, 
by the continuum solvent activation model (Fig. 1); (2) the free energy 

of activation for redox reactions in solution, predicted by the theory, 
bears little resemblance to experimental observations (Fig. 2a, b); 

(3) a piece of evidence in favour of the continuum solvent activation 
theory is that it predicts t'rtfrrate constants for the heterogeneous reac- 
tion from a knowledge of the homogeneous one. This only involves 8 the 
concept of the additivity of the reorganization energy term. This fact is, 
however, not predicted only by the continuum solvent activation model but 
also arises from other models of activation where the interaction between 
the reactants is small. 

There are several reasons for the discrepancy between the con- 
tinuum solvent interaction theory and experiment. (1) In the continuum 
solvent activation theory the change in the energy of the ion in solution 
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leading to the activation required for charge transfer is governed only 
by the interaction of the ion with the solvent, treated as a dielectric 
continuum. This puts out of consideration the inner solvation sphere of 
the ion 6 . Such a model might be true, as a first rough approximation, 
for ions large enough to have no primary solvation sheath attached to them, 
as is clear from Fig. 2a. For reactions involving aquo- complexes, 
however, it is the oriented water molecules in contact with the ion which 
contribute most of the solvation energy of the ion and which therefore 
should make a significant contribution to the energy of activation of the 
ion. Consequently the use of the continuum solvent model to calculate 
energy changes is inappropriate for such reactions. This is also evident 
from Fig. 2b. In such cases, the inner sphere contribution becomes the 
dominant factor. This is supported by the demonstration by Kevan et al 9 
that the energy of a solvated electron is not represented by a continuum 
Born-type expression but principally by nearest neighbour interaction 
with water dipoles. The continuum estimate is in error by about 2 eV. 

(2) It is implicit in the continuum solvent activation model 1 4 
that the energy of the perturbations which gives rise to the activation 
of the ion is proportional to the square of a displacement within the 
solvent. This might indeed be so if the activators were, for example, 
water librators at relatively large distances (for example 50 X) from the 
ion. It is not a likely representation of the energy-distance relation 
if the nearest neighbour interactions of the charge in solution are taken 
into account, because the vibrations and rotations relevant contain 
anharmonicity. It is this lack of anharmonicity in the classical theory 
which is the origin of the poor prediction of the relation between log 
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rate and overpotential. If such anharmonicity is taken into account 
there is excellent agreement between theory and experiment (that is, 
there is linearity between log i and q) for many decades 8 ’ 10 . 

(3) One of the reasons which has been given 6 for neglecting 
the inner sphere contributions in the continuum solvent model was the 
quantal character of the vibrations in the inner sphere. As a recent 
spectroscopic discussion of H^0 + has shown, however, rotational motions 
in solution are classical 11 and there is usually strong coupling between 
vibrational and rotational motion; the energy of the inner sphere may, 
therefore, be activated classically. Together with the classical 
states induced by the Bornian contribution to the total solvation energy 
of ion, the total energy of the initial state will in effect be classical. 
This has an important consequence because no discontinuity in the log i-q 
relations is seen although the structure can be analysed to show up energy 
changes of 0.001 eV. 

One of the remaining uncertainties concerns very large ions: 
does the Landau-Bom expression play an important part as a limiting case - 
that is, if the ion is so large that there is no orienting structure around 
it? For this to be so, it is necessary to assume the applicability of 
the theory of polarons in water; this is, however, unlikely 9 ’ 12 . 

The activating modes in most of the electron transfer reactions 
are therefore not primarily the non- equilibrium polarizations of the solvent 
treated as a continuum. 
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Figure 1. 


Figure 2. 


ANNOTATIONS FOR FIGURES 

3+ 4+ 

Tafel line for redox reactions. (Ti •• Ti ). Theory; 

□ , experiment. 

a. Plot of AG* , against AG* . for electron transfer reac- 
calc ° expt 

tions involving reactants with ligands other than water or 
ammonia; b, plot of A G * alc against A(5* X p t for electron trans- 
fer reactions involving reactants with water or ammbnia 
molecules as ligands. 
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ON QUANTUM ELECTROCHEMICAL KINETICS 
Abstract 

A critical estimation is made of the theory of electrochemical 
kinetics in terms of two models of activation, thermal and electro- 
static. It has been shown that the electrostatic model cannot predict 
the linearity of a Tafel line and also the constancy of symmetry factor 
(3 with overpotential and the relation of rate constants k^ om and k^ 
does not diagnostically support the electrostatic model. 

A quantum mechanical derivation is given to show the con- 
sistency of the time dependent perturbation theory with the GAMOW 
equation for the time dependent barrier concerned in charge transfer. 

1. HISTORICAL 

Few examples 1 ’ 2 ’ 3 are available where an attempt has been made 
to investigate the quantum mechanical aspects of reaction rates, even 

in the gas phase. However, quantum mechanical aspects of electron 
transfer reaction in solution were first discussed by GURNEY 4 in 1931: 

it was one of the first topics in Chemistry to receive quantal discussion. 

The Gurney 4 theory was developed by Butler 5 , Christov 6 , 
Gerischer 7 and Bockris et al 8 . 

The basic aspect of this type of model is that there ^s a ther- 
mal equilibrium between the solvated ions and the solution. A Boltzmann 
distribution is assumed to give a population of excited vibrational states 

of the ion, and the rate of the electron transfer reaction in the product 

of the probability of a quantum mechanical transition (usually of an 
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electron) and that of the probability of finding a suitable acceptor 
(or donor) level in the ion. This model will be referred to as the 
thermal activation model . 

Another model originates with Libby 9 , who suggested that, in 
a redox reaction, the energy change during electron transfer could be 
expressed as the difference between the Bom solvation energies of the 
initial and final states (but see ref. 10). Weiss (1954) 11 , followed 
by Hush (1968) 12 and Marcus (1959) 13 developed a view in which the 
optical polarisation associated with a fast electron transition was the 
principal part of an expression for the heat of activation. In 1965 
Marcus 14 developed a more detailed model, which took into account 
contributions to activation from the inner solvation sheath. The 
former approach, - that based on the Born equation, - has continued to 
be developed by Levich 15 , with Dogonadze 16 and Kuznetsov 17 , with 
explicit stress upon the validity of the application of the Born equa- 
tion and the rejection of contributions by the inner shell of activation. 
This latter model (that which stresses the Born equation) will be referred 
to here as the continuum electrostatic approach . 

Until 1970, the physical meaning of the continuum electrostatic 
view was not differentiated from that of Weiss 11 , Hush 12 and Marcus 13 . 

It did not contain a molecular-level description of how energy was trans- 
ferred from the solvent to the entity to receive charge (which must exist 
in an activated state at the moment of transition). In 1970, Levich 18 

suggested a molecular model which sought to rationalise the continuum 
electrostatic view in terms of fluctuation theory. 
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In this report, we would like to make a comparison of the 
two models*. 


II. THE TWO MODELS 
1. Description of the Models 

The following is common to both models: 

(a) A radiationless electron transfer is assumed. Therefore, 
one calculates a condition for setting up equal electronic energy states 
on the electrode and the solution side of the barrier, so that electron 
transfer can occur from the state Ep, the Fermi level in the electrode, 

to a state of the same energy in the solution. 

(b) As with the consideration of the quantum mechanical 
transitions in reactions, the Born-Oppenheimer approximation is invoked 
to separate the fast moving electronic coordinate and the coordinates 
of the heavy particle, e.g., at ^0* ion, or the solvent. 

(c) In considering the hydrogen evolution reaction, both 
approaches consider the quantum mechanical aspect of a proton leaking 
through the barrier. But in the continuum electrostatic approach, the 
transfer of a heavy particle is pictured as occurring by means of a 
"quantum mechanical transition" of the whole system . In the develop- 
ment of the thermal view the transfer of heavier particles over the 
barrier is taken into account. 

* Importance of the field gains because many biophysical reactions 
probably involve interphasial charge transfer 19 . 
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We now describe the differences. 

(i) Thermal Model 

The model assumes that the ion is in thermal equilibrium 

with the surrounding solvent and there is a Boltzman distribution of 
electronic energy levels in the ion and its solvent shell which provides 
a distribution of levels to correspond to the varying electronic level 
which occurs as the potential of the electrode is changed. Therefore, 
there must be some states in the H^0 + ion which possess an energy to match 
the Fermi level of electrons in the electrode. It is to these states 
that electron transfer from the metal occurs without gain or loss in energy ^ 

In this approach, the permeability of the barrier to protons is 
taken into account. Depending upon the energy level in the proton, at 
penetration, there is a probability of electrons being in the same level 
in the metal, whereupon a transition of the electron to the penetrated 
protons occur. Some protons will also go over the barrier. For heavy 
ion transfer, all the transitions will be over the barrier. For these, 
the quantal character of the electrode reaction will be (effectively) 
restricted to the electronic transition. 

An electrostatic c ontribution to the energy of the initial state 
(ion in the double layer and electron in the metal) was made explicit by 
Parsons and Bockris 20 and Conway and Bockris 21 . They took into account 
the electrostatic energy of the ion in defining the electron energy level 
in the initial state. These energy contributions were divided into chemical, 
ion-dipole and Born terms. They were thought of as contributions to the 
energy of the system at equilibrium. Vibrations of the inner solvation 
sphere will, however, change the radius term in the Born equation and can 
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produce a continuum electrostatic fluctuation, in addition to whatever 
other fluctuations may arise in the electronic energy of the ion due to 
the librations of nearest neighbour dipoles. 


(ii) Electrostatic Continuum Model 

This model differs radically from the basic model used in the 
theory of activation in kinetics, where the small number of excited bonds 
are the medium of the reaction. In the electrostatic continuum model, for 
a bond for which hv ^ kT, the occupancy of the level (in the thermal 
distribution) is too small at room temperature that it does not play a 


part in the reaction. But species have to be excited. To achieve this, 
the reacting ion is pictured as existing in a dielectric continuum, and 
dielectric polarization fluctuations in the continuum are suggested as the 
cause of fluctuations in the energy of the reacting species. 

Levich 18 suggests that the fluctuating particles are water 


molecules, i.e., dipoles, undergoing librations in solution. They have 


a frequency of about lO^ sec *, so that for these ^^1, and they behave 


classically at room temperature* . These librators produce longitudinal 
polarization waves and give rise to electrostatic fluctuations in the energy 
of the ion. Since a large number of librators are needed for the activation** 


-11 


* 10 ** sec is the Debye relaxation time in water. This has been attributed 
to water dipole librations by Levich 1 5i 16 ’ 1 8 . 

** Levich 18 calculates the approximate number in the following manner: If E 


be the activation energy, and each dipole has the capacity of contribut- 
ing an amount of energy not exceeding its own energy hu) Q , then if N dipole 
takes part in creating the reaction situation, the value of N is given as: 

N > E /hu 
, a o 

For E ~ 0.5 ev, the N 10 . 
a 11 

However, using a) = 10 , the frequency quoted by Levich for the librators 


m water, 


hW o - YqTT kT. 

Thus, the majority of librators would possess about 100 quanta, in which case 
the number of particles needed for activities would be about 10. Whether a 




59 


possess, it is necessary that these longitudinal waves are in phase, 
i.e., the oscillators librate in phase to create a standing wave, for, 
otherwise, sufficient energy cannot be transferred to the ion. Thus, 
the concept of the polaron theory developed for solid crystals is used 
to describe the activation process. In this way (though not explicitly), 
it is suggested that energy pulses of the order of 1 ev. reach the ion, 
placing it thus momentarily at an energy level at which radiationless 
transfer can occur. 

Second-order time-dependent perturbation theory is used to 
solve for the quantum mechanical transition probability of the electron 
and for the proton in the hydrogen evolution reaction. 


2. The Rate-Overpotential relation according to both Models* 

Both the models start off from the fundamental equation, which 
gives the rate of an electron transfer reaction as**: 


i 


eC 

s 


n (e f ) P (e f ) W if (e f ) de £ dx. 


( 1 ) 


' S' x 

where i is the current per sq, cm. and gives in fact the rate of the 
electron transfer reaction. 

The main difference in approach is in the calculation of the 
transition probability, W ££ (e £ ). We will examine the methodology in a 

* For list of symbols, see end. 

** The double integral is annulled in the following simple treatment by 
taking the electrons only from the Fermi level and the distance of the 
acceptor states at the outer Helmholtz plane. 
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simplified way to show, the essence of each model . 


(a) The Thermal Approach 

One may consider that electron transfer reactions are adiabatic, 

i.e., the, electronic motions are faster than the motions of the heavy 

particles (the ions) . The probabilities of electron transfer (P^) , 

and of attaining suitable acceptor states (P„ n +(e)), are hence indepen- 

H 2 U 

dent. The transition probability (taken here, for simplicity, as 
occurring only from the Fermx level) is: 


W if C£ P = P T P . 


h 2 o 


C2) 


The problem is to evaluate these two probabilities separately. For the 
hydrogen evolution reaction, the rate -determining step may be chosen for 
the sake of discussion as that corresponding to: 

(H 3 0+) dl + e"(M) + (M - H ... 0H 2 ) dl . 

The potential energy profile will be as shown in Fig. 1. 

The vertical transition (see Fig. 1) AB corresponds to the process 

ef (M) + (H 3 0 + ) dl (M - H ... 0H 2 ) d] _. 

i.e., of taking an electron from the Fermi level of the metal to an 
ion in its ground rotation-vibration state, with no change in the proton- 
coordinate, d. The energy change of the system for this vertical process 
is, say, o AH Q (e), (Fig. 1). Such a transition for an H^O ion not in its 
ground state will be accompanied by an energy change AH Q (e). Eventually, 
at the intersection point X of curves "a” and "b" in Fig. 1, AH Q (e) = 0 
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and electron tunneling from the metal to the proton becomes possible*. 

Since the thermal model assumes that a Boltzman distribution 
exists between the various vibrational-rotational energy states of the 
H + - 0 H 2 ion, the probability of finding the H + - 0 H 2 ion sufficiently 
activated, so that its energy will correspond to the intersection 
point X of Fig. 1, is given by**: 


P + (e) = exp 


h 3 ° 


Ae 
kT ’ 


(3) 


From Fig. 1, Ae is a fraction of o AH Q (e), 


where 0 < g < 1. 

Hence : 


Ae = g AH (e) 
o o 


P + (e) = exp 


H 3° 


g o AH o Ce) 


kT 


(4) 


Since: 


o^o 1 & 


A<}> 


= AH (e) + eAd> 
o o T 


we have 


P + (e) = exp 


H 3° 


B^Ce) 


kT 


exp 


geM 


kT 


(5) 


* For the h.e.r., the r.d.s. of which is (H^O )^ + e (M)-*-(M-H ... 0^)^ 

the point X in Fig. 1 is the point where the equality I + L- <|> = R + A 
is satisfied (see Bockris and Matthews 8 ). Rewriting the above equality 
in the form $=I+L-R-A shows that at the point X the radiation- 
less electron tunneling condition is also satisfied (see Gurney 4 ) . 

** The presentation is heuristic. The quantum mechanical properties of 
the proton are suppressed, for the sake of simplicity in the presentatior 
of principles. Such properties have been taken into account 6 ' 5 . 
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The tunneling probability, assuming a square barrier* is: 

1/21 


P T = A exp 


^{2m(E x - Ep) } 


( 6 ) 


Hence, the total current, using equation (6), is**: 


1 A<f) 


2eC 

8 h 3 


(- ^^-{2 m (Ex - Ep)} j x exp 


3(AH o (e) + eA(j) e ) 
kf 


exp - 


M 

kT 


i o exp 


gen 

kT 


( 7 ) 


This is Tafel's lav;. 


(b) The Electrostatic Approach 

As mentioned in Sec. IV, the electrostatic approach makes an 
initial assumption in contrast to theories of activation in collisional 
kinetics that bonds for which hV/kT ^ 1 do not take part in the reaction. 
The activation is effected by classical degrees of freedom. Moreover, the 
ion-solvent interation in the inner solvation sheath gives rise to the 


* Again, the barrier has been taken here as the simplest for heuristic 
purposes. Christov 6 , and Bockris and Matthews 8b discussed the quantum 
mechanical properties of the proton in terms of Eckert barriers. 

** Ep varies linearly with overpotential. However, E^ also varies linearly 

witlji overpotential. Thus, as the overpotential changes, the state in 
H^O to which an electron is transferred changes by the same amount and 

the height of the electron transfer barrier remains independent of potential. 
There will be a change in field strength in the double layer as a function 
of overpotential. However, the expected changes (say, a doubling of 
field strength) can be shown 80 to have a negligible effect on the shape 
and height of the electron transfer barrier. 
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inequality hv ^ kT. Thus, the number of particles in the higher energy 
states is few. This is usual in any reaction kinetic situation. In 
the present theory, it is assumed that there are too few of such higher 
energy bonds to be of interest in the process. The suggestion is that 
classical harmonic librators in the bulk of the solvent produce (by 
means of long range electrostatic interactions) fluctuations in the 
electronic energy levels of the ion, to give it its required energy 
for a quantum mechanical transition. 

Thus, the potential energy profile in the initial and final 
state are essentially produced from the intersection of two harmonic 
potential wells, as shown in Fig. 2. 

Corresponding to the heuristic presentation given of the thermal 
model, we regard the transition probability as the product of the electron 
transfer probability and the probability of the acceptor state attaining 
the level where a radiationless electron transfer takes place. Thus: 


W if (£ f) = P T P(Ae) ’ C 8 ) 

where now is the transmission coefficient, or electron transfer 
probability, and P(Ae) is the probability that the acceptor state is in 
its activated configuration ready to accept an electron. To evaluate 
P(Ae), consider the PE profile of Fig. 2. Putting the equation of the 
parabolas in the form 24 *: 

U = k(q - d) 2 ) 


and 

U - Q = kq 2 


(9) 


* The assumption of simple harmonic behaviour arises because the displace- 
ment of each oscillator involved in contributing to the activation of the 

-3 

ion is only about 10 ev. 
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where k is a force constant. For the purpose of simplification, we 
have set the k's for the two parabolas to the same. 

At the intersection point, we can therefore write: 

kd 2 - kdq 1 - Q = 0. (10) 

A reorganizational energy, Eg, can be defined as the energy 
difference between the ground state of the reactant and the point corres- 
ponding to its displacement q to overlap the ground state of the other 
oscillator. Hence, by definition (Fig. 2): 

Eg = kd 2 . (11) 


Therefore, from equation (10): 

Eg + Q + 2kq'd = 0 


or 


q' = - 


E S + Q 


2kd 


Hence , 


2 tE S + 
E act ’ >“>’ - — 


4E, 


( 12 ) 


Since the potential energy profiles are made up of classical oscillations, 
the Boltzman distribution is valid. 

Therefore : 


P(Ae) = exp 


E 


act 

= exp 

kT 

/ 


(Q ♦ Eg) ‘ 


4E s kT 


(13) 


i = eCgp o P T n(e £ )exp 
E 


CQ+ Eg) 




4E s kT 


■■= A exp 


s + j =r + 


4kT 2kT 4EgkT 


O’?)' 


Thus 
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The potential dependent part is lumped in Q, since: 


where J^ n and J fn are the ground state energies of the initial and final 
states. Therefore, we can write: 

Q = Q o + en. 

If, in equation (14), 4EgkT, then: 


i = A exp 



ef| 

2kT 


= A ' exp(- jgj. (15) 

Thus, we obtain the Tafel equation. 

2 

However, to obtain this result Q < 4EgkT. This is not usually 
true. Thus, the electrostatic continuum model does not give a linear 
Tafel line. There are no regions on the overpotential-log current relation 
which should be linear, in contradistinction to experiment. 


3. Differences between the Two Models 

The main difference between the two models is in the mode of 
energy communication between the solvent and the ion. 

In the electrostatic continuum view, vibration-rotation levels 
due to the thermal equilibrium of the reacting species are neglected. 
Electrostatic fluctuations from long range interactions of a large number 
of dipoles (see footnote) are envisaged. The electrostatic contribution 
from librations of the inner solvation shell is neglected. 
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The reason put forward 18 for the neglect of vibrations of 
the inner solvation sheath, is that, since the ion-solvent interaction 
of the ion is large, the thermal movements of the solvent molecules in 
the inner solvation sheath hw/kT > It seems to be suggested that, 

because of the small excited to higher states, the effect of vibrational 
activation of the inner solvent sheath can be neglected (See Section IV) . 

Th.e other main difference between discussions in terms of 
thermal and electrostatic activation is in the treatment of the transi- 
tion probability. In the electrostatic continuum view, the treatment 
has not been explicitly made, but it is suggested that it be made by 
means of time-dependent perturbation theory. In the thermal view, a WKB 
approximation is applied to obtain a one-dimensional probability of 
passage through a barrier (Section IV) . 


III. COMPARISON WITH EXPERIMENT 


1. Relation between kj^ and k^^. 


An argument which has been used to support the continuum 


electrostatic theory 


25j 13 


is that the rate constant for the heterogeneous 


reaction can be predicted from the rate constant for the homogeneous 
reaction. The prediction is: 


^hom 

1/2 

M 

z hom 

V. ^ 


Z het 


(16) 


where k, and k, . are the rate constants for the exchange reaction in 
nofR net 


3+ 2+ 

solution (e.g., Fe + Fe 


Fe^ + + Fe^ + ) and at an electrode (e.g.. 


3+ 2+ 

Fe + e(M) -> Fe ). Let us derive this equation from the most general 


assumptions . 
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3+ 2+ 2+ 3+ 

Consider the homogeneous reaction Fe + Fe ■+ Fe + Fe 


as proceeding along the path: 


(Fe 3+ ) 

equiv . solvent 
config. 


(Fe 2 *) 


equiv. solvent 
config. 


(Fe 3 *) 

non -equiv. 
solvent 
config . 


(Fe 2+ ) 
non -equiv. 
solvent 
config. 


electron transfer 


(Fe 2+ ) 

equiv. solvent 
config . 


(Fe 3+ ) 

equiv. solvent 
config. 


rv 2 +. 

(Fe ) 

non -equiv. 
solvent 
config . 


, n 2+, 

(Fe ) 

non- equiv. 

solvent 

config. 


The free energy of activation may then be written as: 

)-(AF 


AF* 

horn 


(AF + AF 

(Fe " 5 )nesc (Fe ) 


nesc 


3+ + AF 2+ 5 

< Fe >esc < Fe >esc 


(17) 


where (AF 


(Fe ) 


+ AF _ ) is the free energy of the activated 

(Fe ) 

nesc v 'nesc 


state and (AF 


3+, 


(Fe ) 


+ AF _ ) is the free energy of the initial state, 


esc 


(Fe ) 


esc 


3+ 


For the corresponding heterogeneous reaction (Fe + e (M) 


2 + 


Fe ), proceeding according to the path: 


(Fe 3+ ) + e~ (M)-*-. 

f, c 3+ i 
(Fe ) 

electron 

f ?+■ 
(Fe 2 *) 

equiv. solvent 

non- equiv. 


non -equiv. 

config. 

solvent config. 

transfer 

solvent config. / 


2 +, 

(Fe ) 

equiv. solvent 
config . 


The activation free energy for the forward and backward reaction can then 


be written as: 
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AF 


S* *" 2*. >- 4F 


CFe- ) 


nesc 


CFe ) 


nesc 


rn 3 +-. 
(Fe ) 


(18) 


esc 


and 


AF. 


hat ’ I' 4 * St. + AF 2+. > - aF 2+. 

(Fe ) (Fe ) (Fe ) 

v 'nesc v 'nesc v 'esc 


(19) 


where tt(AF + AF _ ) is the free energy of the activated 

2 (Fe 3+ ) (Fe 2 ) 

'nesc 'nesc 


state and AF _ and AF _ are free energies of the initial 

state. 


(Fe ) (Fe“ ’ ) 

'esc 'esc 


The two activated states are in resonance and have equal 
energies. Hence: 


I(AF + AF ) = ~ x 2AF = AF 

2 (Fe 3 ) (Fe 2 ) 2 (Fe 3 ) (Fe 3 ) 

v 'nesc ^ 'nesc 'nesc 'nesc 


Similarly, 


I(AF * 4F ) = AF 

(Fe ) (Fe ) (Fe ) 

v 'nesc k 'nesc nesc 


Hence: 


-> 

AF* = AF - AF 

het (Fe 3 ) (Fe 3 ) 

v 'nesc v 'esc 


iF het * iF 2+. - 4F 2t. 

^ 6 ^nesc e ^esc 


( 20 ) 


Moreover, at the equilibrium potential, and equal concentration of 
reactants and products: 
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AF* = AF* 
het lie v. 


( 21 ) 


Now, we can rewrite equation (17) as: 


AF* = AF » AF + AF- 

hom 3 + -» ,-3+. rn 3+.. 

(Fe ) fFe ) (Fe ) 

v ^nesc ^ •'esc v •'nesc 


- AF 


(Fe 2 *) 


esc 


- AF het * iF het • 


So, at the equilibrium potential. 


( 22 ) 


But: 


and 


AF* = 2AF* . . 
hom het 


k, = z, exp 
hom hom 1 


k het - 2 het exF 


AF * 
hom 


RT 


AF* . 
het 


RT 


(23) 


(24) 


Hence: 


Thus: 


hom 


hom 


= exp 


AF* 

hom 


RT 


= exp 


2AF* „ 
het 


RT 




et 


het 


* N 

^hom 

1/2 

\et| 

Z homj 

V. J 


z het 

y 


No assumptions or models have been used in this deviation. 
Thus, relation (16) does not support particularly the continuum electro 


static view 
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2. Predictions of the Free Energies of Activation from a Continuum 
Electrostatic Theory . 

The free enrgy of activation, AF* a 2 C is given according to the 
continuum electrostatic view 18 by: 


AF* . 
calc 


(E s .t Q)‘ 
4E„ 


(25) 


where 


E S = 


' 1 1 l‘ 

1 

1 ] 

2a, 2a„ r 

£ 

£ 

1 2 

✓ 

opt 

stat 

✓ 


(ne)' 


AF*aic were computed for 52 reactions, including both the electron 
*2+ 3+ 3+ 2+ 

exchange (e.g„, Fe + Fe -* Fe* + Fe ) and usual electron transfer 
2+ 3+ 3+ 2+ 

(e.g., Fe + Ce -* Fc + Ce ) reactions. They were made for various 
ligands as well. 

There are few electron transfer reactions for which the temper- 
ature dependence of the rate has been studied. One can calculate an 

"experimental" value for free energy of activation by equating the rate 

' AF* ' 


exp 


RT 


, where Z^ is the 


constant (experimentally determined to expj 
bi-molecular collision number* of the two uncharged species in solution 
per unit volume, per unit time and is given as 14 : 


*The equation (26) used for Z ^ is valid only for gas kinetic collisions. 

Application to solutions would be a drastic approximation. Let us derive 
a crude Z^ for solutions.. Assume a quasi lattice model for a liquid, 

each reactant A and B being surrounded by N nearest neighbours in a 
hexagonal array. The time for a reactant to jump to a nearest lattice 
site is: _ - 

r AB 2 

tCl AB> = 4 • 
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"8tt1cT1 1|/2 2 

[— J r 


(26) 


Using 10 11 for and the experimentally determined k (the rate constant), 
the experimental AF* can be obtained. 

Gaul 

The is plotted against the experimental AF* value in 

Fig. 3. There is no general correlation between the prediction of the 

electrostatic continuum, view and experiment. A few points lie close to 

the theoretically expected line. An attempt was made to categorize the 

reactions on the basis of the nature of ligands and reconstructed the 

plots as shown in Fig. 4a and 4b. In Fig. 4a redox couples with CN , 

Br , Cl and 0 as ligands give results near to those predicted, whereas 

poor agreement is obtained with ligands like ^0, NHg, dipy, Phen. 

The most likely cases in which the electrostatic continuum 

view would have applicability would be for ions of largest radius. 

There, a set of water molecules oriented around the ion and remaining 

attached to it for times much longer than the ions jump time in diffusion 

"primary solvation", would be unlikely 26 . The continuum viewpoint might 

AF* , 

pale 

be more correct. Hence, a plot of should tend to unity with 

expt 

increasing radius. Such a test is shown in Fig. 5, and shows the theory 
to be inconsistent with this test. 


3. Variation of g with overpotential and the Tafel relationship from 
both the approache s 

The continuum electrostatic theory predicts that for the h.e.r. 
in the normal overvoltage region, the transfer coefficient is given as: 
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3 * 


. j_. , - j: 

— + fin 1 in en 

2 2E s " 2E s ' 


(27) 


Experimentally, it has been found 27 ’ 28 that for h.e.r. the transfer 
coefficient remains constant over more than 1 ev. From equation (27) 
this would be so only if J^, - J° n (heat of reaction at the reversible 
potential) = ep. In Fig. 6 the transfer coefficient 3*, calculated 
according to the electrostatic continuum theory, varies linearly with 
overpotential. Using these values of the transfer coefficients for 
Eg = 2ev, the predicted Tafel line is shown in Fig. 7. There is markedly 
no linear region at variance with the experimental line. 

In the thermal model, 3 is: 

1 1 [1-{(<AH?*> - (1 -3)Fd)/d} 1/2 1 

O [ 1 - { (<AH°*> - A o + 3Fp)> /2 ] 


The variation with potential is shown in Fig. 6. The agreement is good. 

Using these 3 values the Tafel line was evaluated using the 

expression: 


2 1 

Ur H 

4frm e.kT 
e 1 

f+W 

I.Y 

[E - E ] 

3 

h 3 j 

w 

e J 

-w 

1 + exp RT 


x exp [-3 (E - E)/RT]dE . (29) 

The equation was derived by Bockris and Matthews 8 and solved numerically 
using a digital computer. The theoretical and experimental lines are 
shown in Fig. 7. The 3 stays constant over a wide range of potential, 
and the Tafel line obtained agrees with experiment quite well. 
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4. Predictions of the Electrostatic (E) and Thermal (T) Models 


Representation 


Comments 


1. (a) The current density is i = Ae F 
exponentially proportional to 
the overpotential. 

(b) Variation of g with (g> o C 
overpotential . 


2. Magnitude of separation 

factors on different metals. S^ T = 3-20 


T predicts experiment. E 
shows continuous curvature. 
Fig. 7. 

T gives negligible variation 
over 1.5 V for her. E shows 
continuous variation. 

Fig. 6. 

No calculation done on E 
model. T model reproduces 
experiment for high r). 


3. Variation of the separation S H/T = f(r)) 

factor with potential. ' 

4. For reaction which involves (log i) = 


adsorbed intermediates, the 
current is a function of the 
heat of adsorption of one of 
the intermediates involved in 
the rate-determining step. 

5. The rate is a function of 
the solvent dielectric 
constant. 

6. The current density is a | 
function of light. 


7. Plot of AF* vs. 

expt 

AF* , . 
calc 


f (AH) 


(log i) n =f(e) 


T reproduces dS/dp better 
than E. Fig. 8. 

For h.e.r. both models are 
consistent with experiment. 
Fig. 9. 


Predicted well by T model; 
and not by E model. 

Fig. 10. 


(i) =K(v-v Not yet diagnostic 30 . 


AF* vs 
expt 

AF* . 
calc 


No correlation on E theory. 
Fig. 3. 


The thermal activation model gives a better agreement with 
experiment than the continuum electrostatic view, particularly in respect 
to the tests (la), (lb), (5) and (7). However, this may arise from 
approximations at present used in the electrostatic continuum model. 


It is therefore desirable to discuss the likelihood of the models themselves, 
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IV. CRITICAL CONSIDERATION OF BOTH 
APPROACHES 

1. The Activation Mechanism 

The poor current potential relationship obtained from the 

electrostatic approach arises from the parabolic nature of the potential 

energy profile. The small perturbations in libratory oscillations of 

the solvent molecules, far away from the ion, must be harmonic. Since 

these oscillations activate the ion, the potential energy profile must 
be harmonic, too. A source of anharmonicity could be the vibrations 

and librations of the inner sphere dipoles, where there is strong inter- 
action between the ion and the dipole. However, this source of activation 
is specifically rejected in the continuum electrostatic approach 18 . 

The mode of activation envisaged by the electrostatic approach 
involves the fluctuations of about a 1000 dipoles (see Section II) "in 
phase". If the librations were not in phase, the large number of dipoles 
needed could not transfer energy to the ion additively at the same 
instance. Polaron theory in solids, where a system of interacting 
oscillators are coupled to one another, has postulated such fluctuation . 
It does not seem to be unreasonable in a solid. However, in a liquid 
without appreciable periodicity, it is questionable whether the principles 
of polaron theory can be applied at all 31 . Thus, Nemethy and Scherega 32 
showed that in water at any instant the percentage of free water ranges 
between 20^50%. It is doubtful that such a liquid will produce a 
significant probability of oscillators in phase. 

Even if there is a group of some 1000 dipoles libtating in 
phase and contributing their energy additively to the ion, there would 
be another (and another) group of 1000 dipoles librating together but in 
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a different phase. The effects would cancel. 

Such a model might be modified to be applicable where there 
is a highly oriented secondary solvation sheath extending to about 20 % 
from the ion. In such a case, librations of the dipoles in the 
secondary solvation sheath might be in phase, and if it extends to 
about 20 % from the ion, the possibility of 1000 dipoles taking part is 
not unreasonable. However, a recent calculation by Bockris and 
Saluja 26 shows that entropy and heat calculations are difficult to 
bring into consistency with the experimental values if the water 
structure is broken past two layers from the ion. Amis 33 has 
discussed a large structure broken secondary layer. The evidence for 
it, however, is ambiguous for it rests upon the values obtained for 
"solvation numbers" from various methods, which are of such ill defined 
significance. 

Coming to thermal contributions to activation from the inner 
layer, the reason for their neglect (see Section II) was that the 
vibrational levels were far apart (hun »kT) , so that there would be 
discontinuity in the Tafel line. 

However, this is at variance with knowledge of the energy 
levels associated with, for example, aq. H,j0 + and, by implication, other 
solvated ions. According to Falk and Giguere 23 , and Falk 34 for liquid 
^0, there is a continuum of energy states from 60 to 3444 cm , Similar 
conclusions follow from 0'Ferrall et al 22 . 

Thus, classical (continuous) modes of heat transfer are 

available in water and, in particular, aq. H^0 + . The appreciable partial 

molar heat capacity (C_) at room temperature for H T 0 + , indicates a broad 

• P a 


range of frequencies giving classical C contributions and provides 

p 

compelling evidence for the above conclusions. 

By analogy to the gas phase, where translational energy is 
converted into vibrational energy (after collisions), it is reasonable 
to assume that the same process occurs in liquid water. The translation- 
al energy of free waters obeys a Boltzman distribution of energy and is 
in equilibrium with the H^0 + ion in solution. The vibration- 
rotational energy levels of the h^0 + ion will also have a Boltzman 
distribution (at least up to a certain energy level) . There will be 
sufficient activation of 0-H bonds to sustain the reaction over the 
rates observed on various metals*. 

The above arguments presuppose that H^0 + ion does not upset 
the equilibrium distribution among the levels by a permeation of the 
barrier to the electrode. This might be the case where the barrier 
thickness for proton transfer is assumed 18 to be 0.5&. It is difficult 
to accept such a thickness, if one is to accept a model of hydrated ions 
in the double layer in which the double layer capacitance is independent 
* For an electrochemical reaction: 

_ kT _ -AH*/RT . 
nF h C S e 

-11 2 

Taking AH* =15 Kcal/mole and i = 10 amp/ cm , C c comes out to be 'v 

-3 2 0 -3 b 2 

10 molecules/cm . Thus, we need 10 molecules/cm in the activated 

state for the reaction to go. Assuming a Boltzman distribution, the 

number of particles in higher vibrational levels is given as: 

XI M -nhv/kT 
N = N e 
o 

for n = 1, N comes out to be 10 4 molecules/cm , and for n = 2 the value 
-2 2 

comes out to be 10 molecules/cm . Thus, we have plenty of molecules in 
Vibrational states to sustain the above rate. 

13 2 

N q is taken to 10 molecules/cm as usual. 


.1 






J. 


1. - -i J 
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of ionic radius. On this basis, the minimum width of the barrier 
for proton transitions is ^ 2.8 X. 

A compelling piece of evidence that an oriented primary 
solvated ion solvent plays a role in the activation process for electron 
transfer reactions, comes from the study of hydrated electrons. The 
absorption spectrum of the hydrated electron shows a strong adsorption 
in the uv region, which is attributed to the ls-2p transition. 

Jortner 36 calculated hv g using the Landau Hamiltonian (that used in 
the continuum electrostatic model), and Fueki et al 37 considered a 
layer of oriented water dipoles around the electron plus the Landau 
contribution. The results of these two calculations, as well as the 
experimental value, is summarized in Table 1. 

This calculation does show that the inner sphere changes 
between the Is and 2p states of the hydrated electron is large. 

There must, therefore, be some rearrangement of the inner sphere (thermal 
or electrostatic) before the electron transfer can occur. 

TABLE I 

Model (AH) . . hv . 

— 'sol. calc. calc 

Continuum 0.91 ev 0.93 ev 

Semi-Continuum 1.95 ev 1.86 ev 

Moreover, if the inner spheres of the reactant and product are 

different, there must be inner sphere activation to satisfy the Franck- > 
Condon restriction of the electron transfer process 9 . Thus, the reactant 
and product activated state, must be in resonance and identical. If the 
inner sphere is different, activation must occur in it before the right 
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1.7 ev 


1.73 ev 
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activated state can be obtained. The explicit omission of the inner 
sphere 18 implies that the inner spheres for reactant and product are 
the same. Applied to H^0 + + e (M) -*• H a( j s + H ^ 0, such a concept is 
unreasonable, and quantum mechanical treatment based on such a model is 
so far from reality that its consideration loses interest. 

The applicability of employing the Born-Landau Hamiltonian 
for treatment of electrons in a polar medium is unlikely 36 , apart from 
considerations made above. The Bom-Landau Hamiltonian arises when the 
electronic adiabatic approximation is used to treat the problem of an 
electron trapped in a polar solvent. This approximation can be applied 
if the binding energy of the trapped electron is much smaller than the 
binding energy of the medium electrons. This is a reasonable 
approximation in a polar crystal where the electron binding energy is 
low. With polar liquids, the situation is different. The binding 
energies are higher, l-2ev. This energy is comparable to the binding 
energy of the medium electrons (4-5 ev); the electronic adiabatic 
approximation is not valid. To calculate the binding energy of 
electrons in polar crystals, Jortner 36 suggests that the additional 
electron and the medium electron should be treated as an equal basis and 
uses the independent particle treatment better known as the self-consistent 
field scheme. The basic difference between these results and those obtained 
with the electronic adiabatic approximation is due to the fact that in the 
self-consistent field scheme the electronic polarisation does contribute 

to the binding energy of the electron. Thus, the potential then becomes 

. 2 , 


of the form 


2r 
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2 
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,-l 

1 
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involved in the 


expression derived from the electronic adiabatic approximation. Thus, the 


reorganisational energy would not be given by the Born-Landau equation. 

It would, correspondingly, not be given by the time-average Born energy, 
but by an equation for the change in the orientation energy on ejection 
of a charge in a polar dielectric. 

2. The Quantum Mechanical Transition Probability 

There are two ways whereby one can calculate the transition 

probability of electrons or protons between two states of equal 

energies 40 , (a) from Perturbation theory; (b) from Tunneling theory. 

Most authors 33 ’ 34 have treated the transion via tunneling theory. 

In respect to electrochemical kinetics, Christov 6 and Bockris and 

Matthews 8 , following the work of Bell 38 and Johnston 39 , have made numerical 

calculations of the transition probability using tunneling theory and 

barriers of various shape. Levich 18 has criticized this approach, not 

only with respect to the electrochemical situation, but in all 
chemical reactions. The barrier created by solvent fluctuations 

is time-dependent. At time t = 0 the system is in its initial state 

and that part of the barrier connect ed with the presence of an adsorbed 

hydrogen does not exist. Similarly, at t = t the system is in its final 

state, and the first part of the barrier does not exist. It is only the 

finite time interval (t) that the barrier forms and disappears as the 

transition occurs. The tunneling theory assumes that the barrier is 

fixed and independent of time and does not change during the course of 

transition. This latter approximation puts into doubt the use of the 

WKB approximation and the Gamow equation in calculating tunneling 

probability. 
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Consider the situation in Fig. 11. When we apply time- 
dependent perturbation theory, we assume initially the electron is 
in region (I) and after an interval of time, due to the action of a 
time-dependent perturbation (which results in the formation of the 
square barrier), the electron finds itself in region (III). There is 
no obvious way of introducing a small perturbation term in the 
Hamiltonian. However, for treating such problems, Bardeen 41 has 
suggested a different approach. Instead of looking for exact solutions 
of an approximate Hamiltonian, we introduce approximate solutions of the 
exatt Hamiltonian. Therefore, we select the following wave functions 
for the electron inside the barrier (Fig. 11). 


T.(x) = a 2 . e 


-k^x 




V 10 = a 2f e 


+k ? x 


(3o) 


Here, T^(x) has to be matched to the correct solution for 

x ^ x 2 , and will decay in the region x x^. On the other hand, 

^(x) has to be matched to the correct solution for x ^x^and decay 
in the region x ^ x^- 

Initially, the electron is in the state represented by 
4' i (x) and we have to compute the transition probability for the electron 
to go into the state ^(x). Let us form a time-dependent solution as 
a linear combination of T^(x) and T^(x) by the usual method. We write: 


-iE.t -iE f t 

f = C(t)¥ i (x)e 1 + d(t)T f (x)e * 


(31) 


Substitution of equation (31) in the time -dependent Schrodinger equation 
gives us finally 42 the effective matrix element for tunneling as: 
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T., = n(H - E.)¥.dx . 
x£ £ 11 


( 32 ) 


This, after manipulation, can be expressed in the form: 

l T if| 2 = -4 i l b 5fi 2 | a «| 2 • (33) 


2 1 2f ‘ 1 2i' 


m 


In Fig. 11 the wave function in the regions I and III can be 
represented as: 


ik . x 

’ r i ’ a u e : + b if e 


-k^x 


ik x -ik x 

^III = a 3f e + b 3f e 


One can therefore obtain b„, in terms of a_, and a_. in terms 

zr 3t 2i 

of a 2i by solving the standard matching problem 43 at Xj and x 2 of Fig. 11. 
Introducing these results in equation (33) , one obtains : 


.4 2 nAV 2,2. 1 2, ,2 " 2lC 2 W 

2 _ h k 2 16k 1 k 3 |a li | |a 3£ | e 

if ' 2 * ,, 2 2. 2 2. 

m (kj. + > c 2 )(k 3 + < 2 ) 


(34) 


Using the Fermi golden rule of second-order perturbation theory, the 
transmitted current can be written as: 


-j - 2lT| y | ^ j. n _ 
3 f h 1 if 1 dE, 


(35) 


dn • 

where is the density of states without spin in the transmitted wave 


and is given as 43 : 


1 2 dn 


m 


3f ' dE, ” , . 2. 5 

f 2iTh k 3 


(36) 


(where |a 3 ^| is the coefficient of the wave function for the transmitted 
wave) . 


The incident current is: 
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Thus, the transmission coefficient is: 


, _ jf _ lbk^k^ _2K 2 a) 

T Ji k l k 3 (k l + K 2 Kk 3 + K 2 D 6 


( 37 ) 


(38) 


where w = x 2 - x^ = barrier thickness, and 

7 1/2 

<2 = (2m(E x - E p )/li } 

Equation (38) has the same form as Gamow's equation. We have, 
therefore, been able to prove that the time -dependent perturbation theory 
treatment of a square barrier penetration problem does give the same form 
of solution as the WKB approximation. (It would break down only in a 
non-adiabatic process) . 

Another argument which might invalidate the use of Gamow's 
equation, apart from that involving the change of the barrier during a 

transition, involves the fact that the top of the barrier fluctuates with 

time. Only if the tunneling time is less than the fluctuation time of 

the barrier, will the usual tunneling expression be valid. Consider the 

situation in Fig. 12. Due to the permeability of the barrier, there is 

a splitting in energy levels in the two wells. Let this splitting be 

<SE, Let and be the wave function in the wells II and IV. Then, 

we can say that, has the energy E q - 6E and - Yj V has the 

energy E q + 6E. 
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Hence, using the time -dependent Schrodinger equation, we can 
solve that find: 

V = expI-iE^/hHO^jj + y^expOiSEt/h] + (Vjj - T IV )exp[-iEt/h] } 
or 

T = 2 exp[-iE Q t/h] {cos^^- 'F jj + isin^j^- Tjy}. (39) 

We have chosen the phases such that the electron is in well (II) 

at t = 0. At t - hir/SE, the electron will be in well (IV). We define 

-14 

t = hm/SE as the tunneling time. If <5E = 0.1 ev, then t = 10 sec. 

If the fluctuation time of the barrier is more than the 

tunneling time, Gamow's equation is unaffected by this cause. For 

electron transfer reactions in aqueous solutions, the barrier is made up 

of vibrational or librational modes of the initial and final states. 

Since the vibrational and librational modes have a frequency about 

10*^sec * and 10** sec * respectively, the fluctuation time of the 

-13 

barrier should be greater than 10 . Thus, the tunneling time is smaller 

than the fluctuation time, for SE = 0.1 ev. 

We have arbitrarily chosen the value of SE to be 0.1 ev. 
Physically, the value of SE is directly proportional to the interaction 
between the electron and the medium. Fran strong interaction, SE is 
generally of the order of 0.TU).3 ev 44 ’ 45 . In our situation, the 
interaction between the electron and the polar media is strong, an 
estimate of 0.1 ev for SE is reasonable. 

Thus, the solution of the barrier leakage problem, using 
perturbation theory is the same as the one obtained by WKB approximation, 
and Gamow's equation should give an accurate estimate of the electron 
tunneling probability. 
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V, SUMMARY . 

(1) Two models exist for activation in electrochemical 
reactions. In the first (thermal , T) there is considered to be an 
equilibrium of vibrational-rotational levels with the translational 
energy of the solvent. These levels are the main source of acceptor 
and donor states for electrons. The other model (the continuum 
electrostatic, E) discards the part played by such levels. Instead, 
the means of energy transfer from the solvent is a long-range electro- 
static fluctuation which is to occur as a result of additivity of 
interaction between a large number of solvent molecules far from the 
ion. 

(2) The rate-overpotential relation is deduced in a heuristic 
way on the basis of the two models. E does not predict a Tafel-like 
law. 

(3) E neglects the electrostatic interaction between the 
librators in the inner solvent shell. However, the ground state of 
these molecules involved librations which would seem to give rise to 
electrostatically origined fluctuations in the ion's energy. 

(4) The relation of k hom and k hgt does not diagnostically 
support E, as previously claimed. 

(5) The calculated free energy of activation at E is not 
parallel to the experimentally observed one in respect to extensive 
solution data (52 systems) . 

(6) T does give a reasonable account of the current -potential 
relation, both in the constancy of g over 1 volt and its reduction towards 
zero at sufficiently large overpotentials. 
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(7) E does not reproduce the effect of change over solvent 

dielectric. 

(8) The transfer of energy from the far out librators to 
the ion depends upon the applicability of polaron theory to solutions 
in water. However, the polaron theory loses its validity in 
structures in which there is no periodicity of the lattice. 

(9) The assumption that there will be a large gap in energy 

between vibrational states in receptor ions, hence a discontinuous 
Tafel line, on model T, is not acceptable. In solution, there is a 

continuum of energy states in, for example, H^O* over a large energy 

range. for water is large at room temperatures: it shows the 

presence of classical modes. 

(10) ls-2p transitions in hydrated electrons cannot be 
rationalised unless the energy of nearest neighbour water molecules 
(i.e., not a continuum dielectric) is accounted for. 

(11) The Born-Landau equation, - the basis of E-, is applicable 
to polar crystals, but unlikely to be valid for electric charge transfer 

in liquids. 

(12) Time-dependent perturbation theory is consistent with 
the Gamow equation for time-dependent barriers concerned in charge- 
transfer. 

(13) Electrostatic fluctuations contribute to the energy of 
ions undergoing transitions with electrons at interfaces. However, the 
origin of the energy of these fluctuations is in thermal equilibrium of 
the surrounding liquid into the nearest neighbour dipoles of the ions. 

Such an energy would add to the states in the electronic-vibrational 
and rotational levels of the ion, in addition to those which are present 



due to equilibrium between solvent and ion. 

(14) Finally, the essential situation is a distinction 

between a kinetics in which collisional activation is the model and 

statistical mechanics gives the distribution of energy and quantum 

mechanics the probability of transfer; versus a mode in which the 

reaction is supposed to occur because fluctuations in energy. It is 

noteworthy that a considerable theory of reactions in solution exist, 
in terms of encounters and statistical mechanics. Either an 

entirely incorrect type of approach has been taken to solution 

kinetics, or the continuum electrostatic (fluctuation) theory is 


invalid. 
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GLOSSARY OF SYMBOLS 


e = Electronic charge 

Cg = Surface concentration of reacting species 
n(e^) = Fermi distribution 
p(£f) = Density of states 

Wif(£f) = Transition probability for the whole system 

Ae = Difference between the energy of the activated state and ground state 
g = Transfer coefficient 

A<j> = Potential difference across the electrode solution interface 

1 - length of barrier 

m = mass of the electron' 
e 


V e f 


= barrier height 


n = overpotential 
E = activation energy 

cLCt 

K = transmission coefficient 


khom = homogeneous rate constant 

khet = heterogeneous rate constant 

Z^m = Collision number for homogeneous reactions 
Zhet = Collision number for heterogeneous reactions 

J^ n , = Energy of the final state 

J. = Energy of the initial state 
in 

r^ = radius of ion 1 

T 2 = radius of ion 2 

r = distance of closest approach 

e opt = °P^^ cal dielectric constant 

e ^ ^ = static dielectric constant 
stat 

D 4n = diffusion coefficient 
AB 
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r A8 = nearest neighbour dipole 
a = Morse constant 

x q = Equilibrium separation between the two Morse curves 
o* 

<A1L > = activation energy for the cathodic reaction 
Dj = dissociation energy 


89 


REFERENCES 

1. K.J. Laidler, "Theory of Chemical Reaction Rates", McGraw Hill, 1969. 

2. S.H. Benson, "The Foundations of Chemical Kinetics", McGraw Hill, 

1960. 

3. G. Kondratev, "Theory of Chemical Reactions in Gases", ADdison Wasley, 

1965. 

4. R.H. Gurney, Proc. Roy. Soc., A134 137 (1931). 

5. J.A.V. Butler, Proc. Roy. Soc., A157 , 423 (1936). 

6. (a) St. G. Christov, Z. Elektrochem. 62, 567 (1958). 

(b) St. G. Christov, Electrochim. Acta, £, 306 (1961). 

(c) St. G. Christov, This journal, 16, 169 (1968). 

7. (a) H. Gerischer, Z. Physik. Chem. (Frankfurt) 26, 223 (1960). 

(b) H. Gerischer, Z. Physik. Chem. (Frankfurt) 26 , 325 (1960). 

(c) H. Gerischer, Z. Physik. Chem. (Frankfurt) 27_, 48 (1961). 

8. (a) J.O'M. Bockris and D.B. Matthews, Proc. Roy. Soc., A292 , 479 (1966). 

(b) J.O'M. Bockris and D.B. Matthews, J. Chem. Phys., 44, 298 (1966). 

(c) D.B. Matthews, Austral. J. Chem., 25_, 231 (1972). 

9. W. Libby, J. Phys. Chem., 56, 863 (1952). 

10. J.D. Bernal and R. Fowler, J. Chem. Phys., 1_, 525 (1933). 

11. J. Weiss, Proc. Roy. Soc., A222 , 128 (1954). 

12. (a) N.S. Hush, J. Chem. Phys., £8, 962 (1958). 

(b) N.S. Hush, Z. Elektrochem., 61_, 734 (1957). 

13. (a) R.A. Marcus. Can. J. Chem., 37, 138 (1959). 

(b) R.A. Marcus, J. Phys. Chem., 67_, 853 2889 (1963). 

14. R.A. Marcus, J. Chem. Phys., 43, 679 (1965). 



- 90 - 


15. (a) V.G . Levich and R. Dogonadze, Dokl. Akad Nauk SSSR 124 , 123 

(1959). 

(b) V.G. Levich, R. Dogonadze and A. Kuznetsov, Dokl. Akad Nauk SSSR 

179 , 137 (1968). 

(c) V.G. Levich in "Advances in Electrochemistry and Electrochemical 

Engineering" Vol. IV, ed. by P. Delahay, 1965. 

16. (a) R. Dogonadze, A. Kuznetsov and A. Chernenko, Usp Khimi 3£, 1779 

(1965). 

(b) R. Dogonadze, A. Kuznetsov and V.G. Levich, Electrokhimiya, 3, 

739 (1967). 

(c) R. Dogonadze, A. Kuznetsov and V.G. Levich, Electrochem. Acta 13, 

1025 (1968). 

17. (a) A. Kuznetsov, Electrochem. Acta, 13, 1293 (1968). 

(b) A. Kuznetsov, Electrokhimiya, 5, 479 (1968) . 

18. V.G. Levich, "Kinetics of Charge Transport", in "Advanced Treatise of 

Physical Chemistry", ed. by H. Eyring, D. Henderson and Y. Jost, 
Vol. IXB . 

19. J.O'M. Bockris, Nature, 224 (No. 5222), 775 (1969). 

20. R. Parsons and J.O'M. Bockris, Trans. Faraday Soc., 47, 914 (1952). 

21. B.E. Conway and J.O'M. Bockris, Proc. Roy. Soc., A248 , 1394 (1958). 

22. R.E. Moore, 0. Ferral, G.W. Koeppl and A.J. Kresge, J. Am. Chem. Soc., 

93, 1 (1971). 

23. M. Falk and P.A. Giguerre, Can. J. Chem., 35_, 1195 (1957). 

24. A.J. Appleby, J.O'M. Bockris, R.K, Sen and B.E. Conway, "The Quantum 

Mechanical Model for Electronic Charge Transfer at Interfaces", 
M.T.P., Volume on Electrochemistry, 1973. 



- 91 


25. N. Sutin, Ann. Rev. of Phys. Chem., 17, 119 (1966). 

26. J.O'M. Bockris and P.P.S. Saluja, J. Phys. Chem., May 1972. 

27. H. Nurnburg, "Studien Mit modernen Technik Z Kinetik schenneller 

chem. and electrochem. Protonen Ubergangs Prozessen", 

Bonn, 1969. 

28. J.O'M. Bockris and A.M. Azzam, Trans. Faraday Soc., 48, 145 (1952). 

29. A.R. Despic, J.O'M. Bockris, J. Chem. Phys., 32_, 389 (1960). 

30. A.M. Braksby and Yu, Ya Gurevich, Electrochem. Acta 1J5, 1245 (1968) 

31. J. Jortner, Conference onMetal Ammonia Solutions, Ithaca, New York, 

1969. 

32. G. Nemethy and H.A. Scherega, J. Chem. Phys., 41, 680 (1964). 

33. J.F. Hinton and E.S. Ami§ Chem. Revs. 71_, 627 (1971). 

34. M. Falk and T.A. Ford, Can. J. Chem., 4£, 1699 (1966). 

35. J.O'M. Bockris, M.A.V. Devanathan and K. Muller, Proc. Roy. Soc., 

A274 , 55 (1963). 

36. J. Jortner, Radiation Res. Suppl. £, 24 (1964). 

37. K. Fueki, D. Feng, L. Kevan and R.E. Christoffersen, J. Phys. Chem. 

75, 2297 (1971). 

38. R.P. Bell, "The Proton in Chemistry", Cornell Univ. Press, Ithaca, 

New York, 1959. 

39. H.S. Johnston, "Gas Phase Reaction Rate Theory", Ronald Press 

Company, 1966. 

40. S.P. Sharma, S.J. Fonash and G.L. Schrenk, Surface Science, 23, 30 

(1970) . 

41. J. Bardeen, Phys. Rev. Lett., 6 _, 57 (1961). 

R.K. Sen and J.O'M. Bockris, in course of publication. 


42. 



- 92 - 


43. H. Eyring, J. Walter and G.E. Kemball, "Quantum Chemistry", John 
Wiley (1944). 




Total potential energy 











kcal/mo le 



cole . 1 



AF^ keal/moie 


Figure 4(b) 


Fexpt. 




0.5 


1.5 


2.0 


1.0 

V (volts) 
Figure 6. 



log i (amp/cm 7 ) 


Figure 7 


-8 









Figure 9(a) 


§ 




10 20 30 40 50 60 


%Methdnot " “» » 




Figure K). 






93 - 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


ANNOTATIONS FOR F IGURES 

Potential energy-distance profile for h.e.r.: is the 

ground state energy level of the reactants; curves (a) 
and (b) show the variation of potential energy with 
intemuclear separation for M(e ) ... H + -OH 2 and M-H. . . 
0H 2 , respectively; q AH o ^ is the standard enthalpy 
of reaction for reactants in taeir ground state. 


Potential energy-distance profile using harmonic oscillators. 


Plot of AF* alc (from electrostatic approach) against AP* X p t 
for electron transfer reactions. 


(a) Plot of AF£ a i c (from electrostatic approach) against 
AF* for reactants with ligands other than water or ammonia. 


(b) Plot of AF* a ^ c (from electrostatic approach) against 
AF*xpt for reactants involving water and ammonia ligands. 


AF* . 
calc 

Plot of against the distance of closest approach for 

expt 

the two reactions . 


Variation of transfer-coefficient with potential. 

Tafel lines from thermal and electrostatic approaches compared 
with experiment. 


Figure 8. 


(a) Variation of the separation factor with potential- 
electrostatic approach . 


(b) Variation of the separation factor with potential- 
thermal approach. 


Figure 9. 
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with the heat of adsorption in 


(b) Variation of 
thermal approach. 
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Figure 10. 
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Figure 11. Tunneling through a square barrier. 


Figure 12. Double well problem. 
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ON THE THEORY OF TUNNELLING IN ELECTRON AND 
PROTON TRANSFER REACTIONS 

Abstract 

The concept of tunnelling in the theory of electron and 
proton transfer reactions has recently been questioned on the ground 
that the situation is a non-stationary one. It has been suggested that 
time-dependent perturbation theory should be applied to obtain the quantum 
mechanical transition probability. We have done this for a square 
barrier. The result for most reactions is the same as obtained by 
the WKB approximation. 

1. INTRODUCTION 

A large amount of work present in the literature (1) concerning 
theories of electron or proton transfer reactions in solution, either at 
an electrode or in the bulk of the solution, uses the concept of electron 
and proton tunnelling through a barrier. For electron transfer reactions 
at an electrode-solution interface, the theory has been developed 
extensively by Gurney (2) , Gerischer (3) , and Bockris and Matthews (4) . 
These authors also consider proton tunnelling at the electrode-solution 
interface. Marcus et al (5) have used the same concepts for electron 
transfer reactions in the bulk of the solution. Bell (6), Johnston (7) 
and Caldin and Kasparian (8) have considered proton tunnelling in the 
theory of proton transfer reaction in the bulk of the solution. However, 
recently (9, 10), the tunnelling approach has been criticized on the ground 
that the process of both electron and proton transfer at an interface 
is of a non-stationary nature. It is contended that the use of Gamow's 
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equation is then erroneous: time-dependent perturbation theory must be 
used. We have calculated the transmission probability of a plane wave 
going through a square barrier which is caused to change by the transition, 
using time-dependent perturbation theory. 

Another aspect of the same problem involves the fact that the 
barrier itself may fluctuate with time, independently of the act of 
transfer, which causes it locally to decompose. Then, also, Gamow's 
equation may be invalid. We have also examined this problem. 

2. CALCULATION 

2.1 The WKB method 

Consider a square barrier as shown in Fig. 1. In the three 
regions the wavefunctions may be represented as: 


^1 = a li ex P + b 1£ expC-ikjX), 

^11 = a 2i ex P( _K 2 X ' ) + b 2f ex P(K 2 x )» (i) 

^III = a 3i ex P( ik 3 X ) + b 3 £ exp(-ik 3 x). 

The well-known solution of the tunnelling problem (11), using the WKB 
approximation, gives the transition probability, P T , as 


J 3 

P T = j ~ 


16k l K 2 


2 2 , 2 2 eXp (" 2K 2 al ^ * 
1 t k i + k 2^ + K 2 } 


( 2 ) 


where w = ^ - x^) is the barrier thickness, and K ^ = { 2 m(V 2 "E)/h } 


2 \ 1/2 
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2.2 Time-dependent perturbation theory approach 

Since P^, is small, a perturbation theory approach seems 
admissible. However, it is not clear that one may introduce a term in 
the Hamiltonian which is small. Bardeen (12) has suggested that, for 
solving problems of this type, instead of introducing states which are 
the exact solution of an approximate Hamiltonian, we should introduce 
approximate solutions of the exact Hamiltonian. We therefore choose the 
following wavefunctions for the electron inside the barrier (fig. 1). 

^ i (x) = a 2i exp(-K 2 x), x<x x 
and 

4> f 00 = a 2f exp(K 2 x), x > x 2 

Hence, ip^ (x) has to be matched to the correct solution for x ^ x 2 » and 
will decay in the region x =< x^. On the other hand, (x) has to be 
matched to the correct solution for x < x^, and decay in the region 
x >x 2 . 

Initially, the electron is in the state represented by ^(x) 
and we have to compute the transition probability for the electron to go 
into the state ^(x). Let us, therefore, express the total wavefunction 
as a linear combination of i|^(x) and ^OO. We write: 

^ I]; = C(t)^ i (x)exp(-iE i t) + d(t)ip f (x)exp(-iE £ t) . (4) 

Substitution of eq. (4) into the time-dependent Schrodinger equation 
gives us: 

C exp(-iE- t)H^ + d exp(-iE^t)H^£ 

= iCi|h exp (-iE-t) + CifK E£exp(-iE.t) 

+ idifh exp (-E.pt) + dijJpEp exp (-iE.pt) . (5) 
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Now, initially, the electron is in ifn(x). Hence, we can 
set C = 1, d = 0 and C = 0. (That C = 0 follows from the normalization 
condition that d(CC* + dd*)/dt = 0). Therefore, (5) gives: 

id^ f exp(-iE f t) = (H - expf-iE.jt), 

or 

id = exp[i(E f - E^t] - E-^-dx. (6) 


Let us, therefore, define the effective matrix element for tunnelling 
to be*: 


if 


^(H - E.^.dx. 


(7) 


We can express the integral in a more symmetric form by subtracting 
iJk(H - E f )i]j*, since this term gives no contribution for the range of 
integration x > x^. 

Therefore, we may rewrite eq. (7) as: 

00 

T if = | ^f (H " E i^i " " E f^f] dx » 


X. < X. < X 

1 B 2 

Integrating eq. (8) by parts gives us: 


( 8 ) 


_ h 

T if " 2m 




dip. 
* __i 
f dx 




dip* 

dx 


= -ihj. 


xB 


if 


(9) 


where is the current operator. Using eq. (3) in eq. (9) yields: 


* This can only be done if the barrier changes slowly with time, for 
example, see Schiff (13) . 
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T if - Ch'K 2 /m)b* £ a 2 . . 


( 10 ) 


One now obtains b 2£ in terms of a 3£ , and a 2£ in terms of a ££ by solving 
the standard matching problem (13) at x.^ and x 2 of fig. 1; one then 
introduces them into eq. (10), and obtains: 


h 4 K 2 16k^kj | a 


2 , . 2 , 


T 1 = 
if 


li 1 


3f 1 


m 


( k i + <9 (kj + 9 


exp(-2< 2 a)) , 


(ID 


where the terms are as defined in eq. (1). Using the Fermi golden 
rule of second order perturbation theory, the transmitted current can 
now be written as: 

"l'~ ,2 ’ (12) 


if = h" x |T if | 2 d n /dE £ , 


where dn/dE £ is the density of states without spin in the transmitted 
wave and is given as (13) : 


dn/dE f = |a 3£ |" 2 m/27Th z k 3 


(13) 


The incident current can be written as: 


3 ± - (hk 1 /m)|a li |‘ 


(14) 


Thus, the transmission coefficient is: 

16k 2 K 2 k 2 
j „ it5K r 2 k 3 


f*T - ^ 2 2 2 2 6xp(-2K_W), 

3 i W k l + 9 < k 3 + ^ 


(15) 


using for j £ the expression resulting from the substitution of eqs. (11) 
and (13) in eq. (12). Comparing eq. (15) with eq. (2), one finds that 
the exponential part is the same, but there is a difference in the pre- 
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exponential part. The ratio of the expressions for P,p given in (2) 

2 

and (15) is k^/k^k^. Since both k^ and k^ have the same order of 
magnitude, the factor will in general be unity. 

Thus, the transmission probability of an electron penetrating 
a square barrier using time-dependent perturbation theory is the same as 
that obtained from the WKB approximation. 

3. EFFECT OF FLUCTUATIONS OF THE BARRIER ON 
LEAKAGE THROUGH THE BARRIER 

An argument which might invalidate the use of Gamow's 
equation both for electron and proton tunnelling, apart from that involving 
the change of the barrier during a transition, involves the fact that the 
barrier may fluctuate with time. If the tunnelling time is more than the 
fluctuation time of thebarr ier, the usual tunnelling expression will 
indeed be valid. To investigate this aspect, consider the situation in 
fig. 2. Due to the permeability of the barrier, there is a splitting 
of the energy levels in the two wells. Let this splitting be 6E. 

Let and be the wavefunctions in the wells II and IV. Then, we 
can say that + ^ IV has the energy E Q - 6E and if) - has the energy 
E q + <5E. 

Hence, using the time-dependent Schrodinger equation, we can 
solve to find: 

ip = exp[-iE Q t/h] + T|>jy3 exp[+i<5Et/h] + - i^) exp[-iEt/h] , 

or 


ij; = 2exp[-iE o t/h]{cos(<SEt/h)4> II + i sin(6Et/h) 


(16) 
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We have chosen the phases such that the particle is in 

well II at t = 0. It is interesting to note that at t = hu/SE, the 

electron will be in well IV. We can hence define t = ^hir/6E as the 

tunnelling time. As we have mentioned before, the tunnelling time t 

has to be small compared to the fluctuation time of the barrier for 

Gamow's equation to be applicable. For electron and proton transfer 

reactions in aqueous solutions, the barrier is made up of vibrational 

or librational modes of the initial and final states. These modes have 

13 -1 

a frequency of about 10 sec for water. Thus the fluctuation time 

-13 

of the barrier for water should be greater than or equal to 10 sec. 

Let us therefore evaluate the tunnelling time. The main problem in 
this calculation is the estimation of 6E, the splitting energy. We will 
use the technique developed by Denison and Uhlenbeck (14) for this 
purpose. These authors found that for a double well potential of the 
form shown in Fig. 2, the splitting energy is given as: 


6E 


h[2m(E - V)l 1/2 

27Tml[exp{(2'irl/h) [2m (E - V)] 1/2 }] 2 


(17) 


where (E - V) is the barrier height, m the mass of the tunnelling particle 
and 1 the barrier thickness. Let us now evaluate the magnitude of <5E 
and the tunnelling time for various electron and proton transfer reactions 
and see when Gamow's equation is affected by the fluctuating character 
of the barrier. 

3.1 Electron tunnelling 

The barrier parameters for electron tunnelling in the hydrogen 
evolution reaction has been estimated for electron transfer from mercury 





4 

■i: 


3 


4 



102 - 


to protons in solution by Matthews (15) to be 7Qkcal/mole for (E - V) 

-8 

and 1 = 2.4 x 10 cm. Thus, <SE from eq. (17) comes out to be 0.15 eV 
and hence 

t = (h/2E)10 14 sec, 

which is significantly smaller than the fluctuation time of the barrier 
-13 

(- 10 sec). Consequently, the use of Gamow’s equation is unaffected 
by fluctuations of the barrier. 

3.2 Proton transfer 

3.2.1 Bulk proton transfer 

Several proton transfer reactions were considered and their SE 
and the tunnelling time calculated. The results are summarized in table 
1. The values of (E - V) and 1 are taken from ref. (8). In five out 
of the six cases, the equation is applicable. 

3.2.2 Electrochemical proton transfer 

From the ground vibrational level of H^O* ion, the height of the 

barrier for proton transfer in the electrochemical hydrogen evolution 

reaction has been estimated by Matthews (15) to be 20 kcal/mole and the 

width to be 2.42 X. However, tunnelling does not occur from the ground 

state of HgO + . The first vibrational level of H^0 + is lOkcal/mole and the 

width is 2 X. The splitting energy for that case using eq. (17) is 0.07 

-14 

eV, and hence the tunnelling time is approximately 1.3 x 10 sec, smaller 
than the fluctuation time for the barrier. Consequently, Gamow’s equation 
is applicable. 
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4, CONCLUSIONS 

In this article we have shown that: 

(i) in respect to the fundamental time-dependence of a barrier 
in a chemical process, the solution for the transmission probability 
according to time-dependent perturbation theory is the same as that 
from the WKB approximation; 

(ii) in respect to the fluctuating effects of the surrounding 
solvent on the barrier, the WKB approximation is applicable in nearly 
all the examined proton and electron tunnelling cases. 
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TABLE I 

DATA RELEVANT TO CALCULATION OF THE SPLITTING ENERGY IN 
PROTON TRANSFER 


Reaction (ref.) 

(E-V) 

I 

6E 

tunnelling 

Comments about 


(kcal/mole) 

(A) 

(eV) 

time (sec) 

WKB approximation 

1. CH 2 Br C • MePh • H+OEt ” 

22.1 

1.59 

0.06 

-14 

1.6 x 10 

valid 

2. RH + D 2 0 a) 

13.1 

1.26 

0.08 

-14 

1.3 x 10 

valid 

3. RH + F" a) 

18.0 

1.17 

0.09 

-14 

1.2 x 10 

valid 

4. RH + F" a) 

24.2 

1.46 

0.01 

-13 

1.1 x 10 

doubtful 

(H 2 0 + NaBr (5M) ) 






5. C 6 H 2 tN0 2 ) 3 C H ; ♦ H0A c 

10.1 

1.66 

0.06 

-14 

1.6 x 10 

valid 

6. H 3 0 + + H 2 0 



0.12 a) b) 

8.47 x 10 -15 

valid 


a) RH = 2 carbethoxycyclo-pentanone. 

b) Estimated by Baker (16) and Somarjai and Hornig (17). 
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ANNOTATIONS FOR FIGURES 

Figure 1. Tunnelling through a square barrier. 
Figure 2. Double well problem. 
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THEORY OF THE PROTON DISCHARGE PROCESS: 

A BEBO CALCULATION 

I. INTRODUCTION 

The theory of proton discharge reaction formulated by Bockris 

and Matthews 1 ’ 2 and Bockris, Srinivasan and Matthews 3 considers the 

activated state to be formed from the excited vibration-rotation levels 

of the ^O-H* bond. The LEP semi -empirical potential energy surface 

had been used to estimate the activation energy, separation factor, Tafel 

line and properties of the activated complex. However, newer criticisms 

of the LEP method by Johnston 4 (difficulty of separating co-ordinates; 

observation of a minimum in the saddle point for H + -> + H) , as 

well as the observation by Christov 5 that the imaginary frequencies depend 

on the coulomb-exchange energy ratio, make an alternative approach to 

proton transfer calculations desirable. 

The BEBO method 6 , on the other hand, is a fully empirical method 

(its empiricism depends on spectroscopic and thermodynamic parameters), 

not open to the criticisms made of the semi-empirical LEP method. This 

method has been developed in detail both for proton and hydrogen atom 

transfer processes in the gas phase with good success 7 . Some attempts 

have been made to use such a method for proton transfer in solution 8 . 

n. THE BEBO METHOD 

The BEBO method assumes that along the reaction co-ordinate the 
sum of the bond orders of the bond being broken and the bond being formed 

is always unit. Thus, for the rate-determining proton discharge step 
of the hydrogen evolution reaction on Hg, which may be written asi 
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H 2 0-H + e(Hg) -*■ H 2 0 


H 


. e (Hg) + H 2 0 + Hg - H ads , 


the energy along the reaction co-ordinate may be written as 11 *: 
E = E° (l-n p ) - E° (l-n)« * BE° [n Y (l-n) Y ] , 


xn 

-,0 „o 


J fin' ' rep 

where ET , E^. , E~ are the energies of the initial (H„0 + + 

in’ fin’ rep 6 v 3 soln 

e metal )j £inal (H ads + H 2 0;) and re P ulsive ( H S 
is the bond order of the final state of index q. 

p and q are defined as 9 : 


( 1 ) 


0) state**; n 
The parameters 


P = 


0.26£n(E^ n /e x ) 

(R - R? ) 
v xy in' 


and 


q = 


0.26£n(E°. n /<^) 

(R* - R° ) 

v xy fin' 


where R , R' and e and e' are the internuclear separation 
xy xy xy xy 

and bond energy respectively in the corresponding noble gas 


diatomic cluster where n = 0. 


R and e are for the initial state 
xy xy 


and R* and e* are for the final state, 
xy xy 

y and B are defined as: 

Y = 0.26a 

and B = exp[-aAR] , 

where a is the anharmonicity constant of the bond Hg-0 and AR = [R? n + R° in 
- R° g p] , where R? n> R ££ n and ^rep are t * ie equilibrium internuclear separa- 
tion of the bonds M-H and Hg-0 respectively. 

* The initial state energy is submitted from the total energy in eqn. (1). 

** By the repulsive state we mean the energy of interaction between Hg 
and the oxygen of the water molecule. 


10 9 


(a.) Estimation of the various parameters of equation (1) 

fi) Estimation of E? , E° and E° : 
v * in fin rep 

E? and E~. were taken from the estimation of Parsons and 
m tin 

Bockris 10 , and Bockris and Matthews 1 . They were -69kcals/mole 
for E? n and -75kcals/mole for E^ n at the reversible potential. E° g p 
is the dissociation of the Hg .... 0 bond, and was taken as 80kcals/mole 
from the thermodynamic estimation of Solomon, Enke and Conway 11 . 

(ii) Estimation of y and B : 

The anharmonicity constant a is estimated from the equation: 


a = (F° / 2E° ) 1/2 

rep rep 


E° is obtained from Badgers Rule. The value of ot arrived at is 
rep 6 

1.49& _1 . Thus, 

B = exp[-aAR g ] = 0.7335 

where R? , R ~ • is estimated from Bockris and Matthews 1 , and R° 

in’ fin * rep 

is taken from Cotton and Wilkinson 12 . 

Moreover, y = 0.26a = 0.3891. 

(iii) Estimation of p and q : 

These are estimated from the equation given in the previous sec- 
tion and the values obtained are p = 1.3and q = 0.8. 


(b) Calculation of the activation energy at the reversible potential and 
the variation of the transfer coefficient with potential 

Since p and q are both very near to unity, we can expand the first 
and second term of equation (1) in a Taylor series expansion, as suggested 
by Marcus 8 . Thus, we get: 
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E=(l-n)(E° n -E° in )-jCp-l)nE° n 5,n n - E° in (q-l)l-n(Sn(l-n) + Be° ep n y (l-n) Y . ( 2 ) 

The activated state corresponds to the maximum of the E vs n curve given in 
equation (2). Thus, we differentiate E w.r.t. n and set it equal to zero. 
Hence: 

g| = -Q-(p-l)E? n (l+An n)+(q-l)E° in [l+£nCl-n)]+yBE° ep n Y ' 1 (l-n) Y ‘ 1 Cl-2n)=0, (3) 


where Q = (E? - E° ). 

m fm J 

Solving the above equation numerically gives us the value of n 
at the activated state. Let us denote this value of n as n*. Thus, 


the activation energy now becomes: 

~>Q n._. _ j. /n j. T-.0 


E* = Q(l-n*)-(p-l)n*E? n &n n*-(l-n*)E^ in (q-l)JlnCl-n*)+BE° ep Cn*) Y (l-n*) Y 

The transfer coefficient 3 is defined as the variation of the activation 
energy with respect to the reaction heat, Q, where Q = Q + eri, and q is 
the overpotential . Thus : 

^{p« (l-n*)-(p-l)n*to n*+ (p-l)E? n (l+J!*i n*)+(q-l)E° in [l+An(l-n*)] 


(4) 


+ ^BE° ep (n*) Y_1 (l-n*) Y-1 (l-2n*)] . 


(5) 


From (3) in (5) : 


dAE* 

dQ 


= (l-n*)-(p-l)n*£n n* = (l+n* p ) . 


( 6 ) 


Since, as the overpotential changes, E? n changes*, the value of n*, AE* and 
3 also change. Thus, using equations (3), (4) and (6), we are able to 
get the activation energy at the reversible potential, and the value of the 


* E? n comprises the energy of the H^0 + ion in the double layer and the 
energy of an electron in the metal, and is therefore potential dependent. 
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transfer coefficients for various overpotential. These results are 
summarized in Table (1) and the variation of 3 with overpotential is 
shown in Fig. 1. 


(c) Calculation of the Quantum Mechanical Rate of the Reaction 

Following Bockris and Matthews 2 , we can write the quantum mechan- 
ical rate of the proton discharge reaction as: 


i = 


oo 

= k l c +• [ 

H„0 J 


W(E)e 


- (E-E°)/kt dE 


(7) 


3 E° 


where k is a constant (not a rate constant) , C is the surface 
1 HO 

+ 2 A 

concentration of H^O ion per cm , W(E) is the proton tunneling probability, 
and E° is the zero point energy of the H 2 0-H + bond. The corresponding 
classical current may be written as: 

i r0 = k i c k T exp[-(E*-E°)/kT] . (8) 

ut i H ^ Q + 

Thus, we can define: 

OO 

r \ r 

exp[-(E-E°)/kT]W(E)dE. (9) 


T = 


1 
1 C£ 


exp (E* - E°)/kT 


Let us define: 


J( l " ) 


W(E) exp [ - (E - E°)/kT]dE . 


( 10 ) 


Eo 


Thus, equation (9) becomes: 



i 



T = 

k J 

j 

" kT 6X P 

(E* --- E°)/kT 

► - ' 4 


V 


cm 


or 
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= x 


cr 


1 

kT 


exp (E* 


E°)/kT v 

. 


( 12 ) 


To evaluate the current one must know the shape of the barrier 
along the one dimensional normal mode. The Eckart barrier appears to 
have the closest fit to the real barrier 2 . It has the form: 


V (x) 


A exp[27tx/d] , B expr27rx/dl 
1 + exp[2ux/d] (1 + exp(2irx/d}^’ 


(13) 


where 

A = 

Q 0 + en , 

(14) 


Qo 

_ E (+oo) " E (-oo) » 

(15) 


B = 

2E* - A + 2{E* (E* - A)} 1/2 , 

(16) 

and 

E* 

= E* - 3e p > 
o o ’ 

(17) 


where 2d is the barrier width and E* is the barrier height at p = 0. 
Using equation (6) in equation (17) : 



E* = E* - (l-n* P )e Q p . 

(18) 

For this barrier: 



W(E) 

Cosh2ir(X+y) - Cosh 2ir(X-y) 
Cosh 2 tt(A+u) + Cosh 2710 ’ 

(19) 

where 

X = (d/h)(2m e££ .E) 1/2 , 

(20) 


y = (d/h)‘ (2m e££ (E - A) } 1/2 , 

(21) 

and 

if,* ,2 B , Jl/2 

° ~ 2 C8m eff d ,2^ ~ 1 

(22) 


h' 


where m e££ is the effective mass of the tunneling particle along the 
reaction co-ordinate. 

For barrier widths greater than 3A and E* > 1.0 x 10 ergs 
we use the approximate formula: 
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wfFi - ex H 2-TrO+u-g) }-exp 2Tr(A-y-cr) -exp{-2TT(A-y+a)} 

^ ' 1 + exp{27T(X+y-a)} ' 


(23) 


(d) Effective Mass of the Tunneling Particle : 

For a particle tunneling along the reaction co-ordinate, the 
system having a fixed centre of mass and constant angular variables, it 
was shown by Johnston that the effective mass of the particle along the 

reaction co-ordinate is different from the actual mass of the particle. 
Thus, Johnson 9 shows that the tunneling particle has an effective mass, 

S 


given as: 


m 


= m A m B (1+c2) ~ m H m B c2 * m H m A 

0-p *P o 

(m^mg+ii^) (1+C j 


(24) 


where 


dr 


C = 


A-H 


dr H-B * 


(25) 


for a system of the form [A...H...B], 

Christov and Georgiev 13 have shown that under the limiting 
conditions , and m H <<m B are °^ ta i nec i from equation (24) . 


m eff 2 m H * 


(26) 


In the system we are considering A is H o 0 and B is Hg and the inequality 
m H <<m A and m H <<m B evidentl X satisfied. Thus, in evaluating W(E) in 
equation (23) , we use equation (26) to define the effective mass of the 
tunneling particle. 


(e) The Separation Factor and its Variation with Potential 

The hydrogen tritium separation factor is defined by the equation: 
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S T ' 


fC H l 


J H 


where and (C^/C,^) g are the ratios of atomic concentration of 

hydrogen of tritium in the gas phase and in solution. The ratio 

(C H /C T )g is equal to twice the ratio of the velocities of H 2 and HT 
evolution since »C^. Using the relation S^, may be expressed 


as: 


t f T f 

S T - a. Ji . -SI . J™* 


where a is a ratio of symmetry numbers 


-T- • f • Kj, , 

T fT r H 2 0,g 


= S, 


T, class x T ’ 


(28) 


Bockris, Srinivasan and Matthews 3 showed that: 


HTO , g^H ? 0 , g^ 316.46. 


It is assumed that the activated complex H 2 0. . . ,H. . . ,Hg (or its isotopic 
analogue) is similar to a linear triatomic molecule (cf. Parsons and 
Bockris 10 ). Thus, ft/f| is given by: 


ft ft ft 

H _ Vh 

ft" ?r" ’ ft ■ pr 

T T,t r,T vib, T 


ft 

vib,H 


(29) 


Since the activated complex is regarded as immobile, the translational 
partition function ratio is unity. Bockris, Srinivasan and Matthews 3 
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estimated the ratio of the rotational (around the axis) partition 
function ratio is: 


f+ 

H,r 


f, 


T,r 


4 


= 0.962 


where the I's are the amounts of inertia of the activated complex. 

Moreover, they wrote the vibrational partition functions as: 


4ib,H Sin h ( H V 2kT ^s Sin ^ 2 ( h V 2kT) b 

Tf ~ = Sin h (hv*/2kT) c * Sin h“ (hv u /2kT). * 
vib , T H s H b 


(30) 


The suffixes s and b stand for stretching and bending frequencies. 
It was further shown that: 


Sin h 2 (hv T /2kT) b m H x 

Sin h 2 (hv H /2kT) b m T 3 


(31) 


For calculating of the stretching vibrational frequencies, it is 
necessary to solve the secular equation: 


X 2 -X 



’1 

h i_] 

F + 

fi 

1 1 
h — — 

F 

- i-F 

. m l + m 2 + m 3 

r 7 

V P 


m l 

m 2 

F n 

sT 

m 3 

22 

m 2 12 

m l m 2 m 3 

F 11 F 22 F 12 







M 


. 

( 


= 0 , 


where is the force constant of the bond between m^ and m^, F^ is 
the force constant of the bond between m 2 and m^, and F 12 is the coupling 
force constant. X is given by the expression: 

X = 4 ttV . (33) 

Bockris, Srinivasan and Matthews 3 used an LEP surface to estimate 
the force constant. We will estimate it from the BEBO calculations. 
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In the BEBO method, if we define the force constant along the 

reaction co-ordinate to be and that perpendicular to it to be F^, 

then F and F are given as 9 : 
p o a 


F p~ 


and 


10.7 


1 + 1 
~~~2 ,, 2 

. (n*) 2 ' p (l-n*) 2 ' q 

n* (1-n*) 


F? n 3 + F° (1-n*) 3 + (F 1 
„ _ in fxn v * K : 

i — 

a 

(n ■' 2 + (n-n* 


E° 2B 


+ Te v.. 


0 / 
rep 




i + (ill) (l-2n*j 
2n*(l-n*) 

(34) 

*) Y B 

, (35) 


where F? , F° and F° are the force constants of the initial, final 
m fin rep 

and the repulsive states. 

To solve equation (32) we have to have the force constants in 
valence bond co-ordinates, with the force constants F^, and F^* 
The following realtion holds 9 between F.^, F^ 2 and F 1<7 , and F^ and F^. 


12 1 


F (1-n*) 2 + F n* 2 
F = -P- g 

11 2 2 

(n*) + ( 1 -n*) 

F (n*) 2 + F (1-n*) 2 

p = p q 

22 2 2 

(n*r + ( 1 -n*)^ 

_ C-Fp + Pp)n* (1-n*) 

^12 ~~ 2 2 
(n*) + (l-n*r 


(36) 


(37) 


(38) 


Thus, knowing n*, F^ and F^ we obtain F^, F^ 2 and F^* Thus, we solve 
equation (32) to obtain the vibrational frequencies of the activated state 
and the ratio of the vibrational partition function given by equation (3) . 
The results are tabulated in Table (2) . 
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(f) Estimation of t^/t^ and the Tafel lines: 

Equation (10) was solved by doing a numerical integration, taking 
the expression given by equation (23) for W(E) on an IBM 1130 digital 
computer for various values of the width of the barrier. The dependence 
of the results upon the width of the barrier is shown in Fig. (2). 

The most probable parameters were selected on the basis that they give 
good agreement with (a) Tafel slope, (b) Activation energy, and (c) 

The variation of the separation factor with potential. Thus, the Tafel 
lines obtained for the most probable barrier parameters are shown in 
Fig. 3a and b and the data are summarized in Table 3. The variation 
of the separation factor with potential is shown in Fig. (4). The 
theoretical calculation was not pursued at higher overpotentials due 
to limitations of the programme. However, the calculations are bound 
to show a constancy of s with q at sufficiently high q's, because when the 
barrier has been made sufficiently low by increasing q the quantal con- 
tribution becomes negligible. 

III. CONCLUSIONS 

Thus, in this article, we have shown that the BEBO method applied 
to the electrolytic hydrogen evolution, gives: 

(a) Excellent agreement on the variation of the transfer 
coefficient with potential. 

(b) Good agreement with the experimental Tafel slopes. 

(c) Tafel lines which show the occurrence of barrierless 
discharge at high overpotentials. 

(d) rGood agreement between the theoretical and experimental 
separation factor values and their variation with potential. 
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TABLE I 


VARIATION OF THE TRANSFER COEFFICIENT WITH POTENTIAL 


0 

(AH*) calc 

C^expt 2 

3 


(kcals/mole) 

(kcals/mole) 


0.0 

15.54 

20 

0.561 

0.25 



0.558 

0.5 



0.554 

0.8 



0.551 

0.9 



0.547 

1.0 



0.591 

1.25 



0.531 

1.5 



0.52 

1.63 



0.48 
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TABLE 2 


„ ^HTO,g 

*T*f 

H 2 0.g 

fT 

r,H 

?“ 

r,T 

F 

dyns/cm 

F P 

dyns/cm 

W 

H 13 

xl0 _l 

sec 

W T 

x!o« 

sec 

^vib.H-’b 

^vib.T^b 

C s ji ) 

i class 

316.46 

0.962 

3.65X10 5 

-2.6xl0 5 

6*35 

1.53 

0.021 

3.6 


120 - 


TABLE 3 


TAFEL LINE AND THE VARIATION OF SEPARATION FACTOR WITH POTENTIAL 



E* = 1.6 x 
o 

10 ergs , 

A = 0.6 x 10" 12 
0' 

ergs, and 2d = - 

4.5X . 


n (in 

ergs) n 

(in volts) 

Jq(H + ) 

Jq(T + ) 

T 

H/T 

S H/T 

0.00 







0.4 x 

IQ' 12 

0.25 

0.5857xl0~ 12 

0.1386xl0 -13 



0.8 x 

10- 12 

0.5 

0. 1080xl0~ 9 

0.3251xl0 -11 

29 

10.15 

1.2 x 

10” 12 

0.8 

0.1843xl0" 7 

0.5442xl0~ 9 

2.4 

8.40 

1.4 x 

10" 12 

0.9 

0.2150xl0 6 

0.6769xl0" 8 

2.25 

7.85 

1.6 x 

10" 12 

1.0 

0.2228xl0" 5 

0.7483xl0" 7 

2.07 

7.24 

1.8 x 

10" 12 

1.13 

0.2455xl0“ 4 

0.8866xl0" 6 



2.0 x 

10" 12 

1.25 

0.2166x10 

0.8399xl0" 5 

1.55 

5.4 

2.4 x 

10- 12 

1.5 

0.1998xl0 _1 

0.1148X10" 2 

1.2 

4.20 

2.6 x 

10" 12 


0. 4091x10 _1 

0.49 80x10" 2 



2.8 x 

io" 12 


0. 4300x10" 1 

0.2769X10" 1 





t ] in volts 
Figure i. 

***- 


M 


1 .8 2.0 




-b (mV) 







m 




-A — At Expt. points, 
-©-©-o- Theoretical points 


(£> 


d'TJ. Theory 


- 6.00 volt 


dl 7 Expt. 


o' 9.00 volt 


Discrepency r 33% 

\ between theory & expt 

p. - 


N 
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Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


ANNOTATIONS FOR FIGURES 
Variation of 3 with overpotential q. 

(a) Tafel lines for various barrier parameters. 

(b) Variation of the separation factor with potential. 

Variation of Jq (the quantum mechanical correction to 
Tafel line), for most probable barrier parameter. 

Variation of separation factor with overpotential for 
most probable barrier parameters. 
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ON ALTERNATIVE ACTIVATION MECHANISMS IN ELECTRON TRANSFER 
REACTIONS IN SOLUTION 

Summary 

The model for electron-transfer kinetics in solution is 
considered. In one model, the appropriate energetic condition 
for charge transfer is met by a small number of vibration-rotation 
states, in thermal equilibrium with the solution. Collisional 
activation (CA) between ion in the solution and the solvent is the 
origin of such states. In another model (EF), CA is neglected, and 
the appropriate energy states are regarded as being reached by the 
fluctuations in the energy of the ion, as a result of its inter- 
action with many surrounding solvent molecules. The methodology 
of deduction of the dependence of the charge-transfer rate upon the 
interfacial potential difference for the two models is outlined. CA 
suggests a linearity of log R with q 0* = rate; q = overpotential). 

There is no linear region suggested on the EF model, although such 
relations are regularly observed. This lack of consistence arises 
because the perturbations which are regarded as the origin of activation 
in EF are simple harmonic. This is an intrinsic part of the model of 
EF, not as an approximation within it. Correspondingly, a comparison 

OX 

of (AF ca j c ) E p with experiment shows inconsistencies. 

Comparison with spectroscopic data for 1^0 H in solution suggests 
that the energy distribution in the vibration-rotation levels in this 
ion is continuous and that classical modes of vibration exist in water. 

A supposed discontinuity, - which would have annulled the deduction 
of Tafel’s law, - was an origin of the EF model. 
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In the EF model, the applicability of the Born-Landau equation, 

F° = e^/2r(l/e - 1/e ) is assumed. However, the applicability 

Op L St 3X. 

depends on a sufficiently large difference of the energy of an electron, 
trapped in the medium and bound to atoms in At, This difference is 

great if the medium is a solid, but not if it is a liquid. 

States suitable for acceptance or donation of electrons from 
ions to metals arise (at the equilibrium potential) much more frequently 
as a result of the equilibrium of the ^0* ion with the solvent heat 
sink than those by electrostatic fluctuation. 

I. INTRODUCTION 

A quantum mechanical formulation of the rate of electron transfer 
reactions in solution was first given by Gurney 1 (1931). This model 
(termed here the collisional activation model (CA)) was further developed 
by Butler 2 , Christov 4 , Gerischer 3 , and by Bockris and Matthews 5 , 

The essence of this approach is that a thermal equilibrium between the 
vibration-rotation levels of the solvated ions and the thermal energy 
of the solvent is assumed. 

Weiss 6 , and Platzman and Franck 7 , founded another model which they 
assumed that the energy of activation arises from some undefined movements 
in the continuum solvent. It was elaborately developed by Marcus , 
considering non-equilibrium dielectric polarization changes in a continuum 
solvent. A molecular model for this type of activation process was 
described in detail by Levich 9 , and by Dogonadze 10 , where it was stated 
that the fluctuations in the electrostatic energy arising from the 
libratory movement of the solvent dipoles cause the activation of the 
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reactant ion. This model will be referred to as the Energy 
Fluctuation (EF) model. 

In this article, we attempt to distinguish between these two 
models. The collisional-activation model is an application to electron 
transfer kinetics of the concepts of thermal distribution of vibrational 
states, as discussed in the reaction rate theory of gas phase reactions. 

The electrostatic fluctuation approach differs fundamentally fxem 
previous models of reactions in solution. 

II. THE TWO MODELS 

In both approaches a radiationless electron transfer is assumed 
and the Bom-Oppenheimer approximation is invoked to separate the co- 
ordinates of the fast moving electron and the slow moving heavy 
particles, e.g., the solvent molecules. The continuum electrostatic 
approach uses a double adiabatic approximation by which it separates 
the co-ordinates of the electron, the proton, and the solvent molecules. 

In considering the hydrogen evolution reaction, both models 
calculate the quantum mechanical probability of protons leaking through 
the barrier. However, they differ in the methodology of evaluating the 
transition probability. In the thermal approach, the WKB tunneling 
theory is used 5 , whereas in the electrostatic continuum approach, the 
use of time -dependent perturbation theory is advocated 9 , although the 
barrier has usually been considered transparent to protons. 

In the CA model it is assumed that there is a classical distribution 
of vibration-rotation levels in the ion-solvent complex (because of the 
increase in the number of such levels in solution, compared to that in the 


126 - 


gas phase). Furthermore, for a reaction like the hydrogen evolution 
reaction, the CA models assumes that the H^O* ion in the excited 
vibration-rotation levels produces the "activated state" for the 
reaction. In the "activated state", quantum levels for electrons 
exist in the H^O* ion which are equal in energy to the energy of 
mobile electrons in the metal, and consequently a radiationless 
electron transfer occurs. 

The quantum properties of the proton (in reactions such as H^O + 
e ** ^ads^ are accounte d f° r by assuming that protons have a probability 
of permeating the energy barrier between it aid the metal. A Gamow 
type lunneling equation is used for the probability of penetration*. 

At present, theories of chemical reaction in the gas phase 
assume that the small number of excited bonds are the medium of the 
reaction. (The type of model is assumed in the CA approach) . 

Exactly analogous concepts are used to explain the mechanism of reac- 
tions in solution. However, the EF model differs from these available 
theories. Here, it is assumed without argument that the number of 
excited levels is too small (at room temperature) to be kinetically 
significant. For a radiationless electron transition, however, there 


* An objection 10 has been made to the use of a one -dimensional WKB 
approximation in calculating a transition probability. Thus, it is 
claimed that the WKB method calculates transition through a stationary 
barrier. However, the transitions involved in the proton discharge 
process are non-stationary in nature. Time-dependent perturbation 
theory should clearly be used. However, it can be shown 17 that, when 
it is applied, then, for the time domain concerned in chemical kinetics, 
it gives approximately the WKB approximation. 
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must be an excitation of the electronic energy level in the acceptor 
so that it becomes equal to the energy of the Fermi level in the metal. 
To achieve this, the electronic energy of the reacting ion is pictured 
as experiencing a fluctuation in its energy with respect to the 
solvent dipoles in the bulk of the solution. 


m . A DEDUCTION OF THE RATE-OVERPOTENTIAL RELATION ACCORDING 

TO BOTH MODELS 

1. The Col lisional -Activation Model 

Considering the discharge reaction, 

H 3 0 + + e - ^ ■> MH + H 2 0, 
the quantum mechanical rate may be written as: 


i 


eC H + 


n(e i )pCe i )W if (e i ,x)de i dx 

x e 


> 


CD 


where C is the surface concentration of the proton, n(E.) is the 
H 1 

Fermi distribution of electrons in the metal, p(s^) is the density of 

states and is the transition probability for the whole reac- 

tion. Now, W- (e_) using the adiabatic approximation, may be split 
into two parts: (a) the probability of electron tunneling, P(e^,x); 
and (b) the probability of proton transfer across or over the 
barrier, W(e^,x). Thus, we can rewrite equation (1) as: 


i = eC 


H 


n(e i )p(e i )P(e i ,x)W(e i ,x)de i .dx . 


X £ 


( 2 ) 
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If we make the assumption that the transition occurs from a fixed 
distance essentially (i.e., the OHP) , we can get rid of the integration 
over x. Hence: 


i = eC 


H 


P(e i )nCe i )pCe i )WCe i )exp 


f- (e i - e o 


)/kT 


de. 


(3) 


where WCe^) is the probability of proton tunneling at an energy level e^, 
e is the zero point energy level for the H~0 + ion. 

O j 

In the adiabatic case, where the electron transfer probability is 
unity, the equation reduces to: 


i = eC 


H 


r 

f(e, - 

n(e i )p(e i )W(e i )exp - 

✓ 

1 0 
kT 


de i * 


C4) 


Using the free electron approximation, for metals, we can write: 

mV e 1/2 

p(e i> ■ -n (2, ( E i - u 0 ’ } ■ 

2 IT h 


(5) 


where is the average volume occupied by one electron, and n(ep we 
substitute the Fermi distribution: 


n(e i ) 


exp 


Ce i-^ +x r 


kT 


( 6 ) 


Using (5) and (6) in (4) we get: 

1-1. mV 


i = eC 


H 


e i ~ e F 

exp (— k ---) + l 


^{2m(e.- e ^Ce^exp 


' 2A 3 ' ^ 0 


(e i _e o )/kT 


de,. (7) 
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A problem is to evaluate W(sz), the proton transfer probability. It 
depends on approximating the actual potential energy barrier by a 
barrier of similar shape for which W(e^) has a known solution. 

The Eckart barrier 5 resembles such a barrier. The unsymmetrical 
Eckart barrier is given by: 


vm = A exp (2 x/d) B exp (2 x/d ) 

1 J 1 + exp (2 x/d) f ...2 ’ 

' 1 + exp (2 x/d) 


( 8 ) 


where 2d is the barrier width, A is the energy difference between the 
potential energy minimum, and B is defined by: 


B = 2e^ - A + 2 


e. . (e. - A) 
x v 1 1 


For the hydrogen evolution reaction (HER) we may write: 


A = A = e 
o o 


£= £ 0 + 3e o n . 


(8a) 


For the barrier represented by (8), the tunneling probability is: 


where 


Wfp > - Cosh 2 tt(A + T-0 - Cosh 2ir(A - y) 
^ z* ~ Cosh 2tt(A + y) + Cosh 2ira 


( 9 ) 


X = (d/h) (2me.)' 


y = (d/h) 


2m (e^ - A) 


a = (h) 




|^(8md B/ 2 ) - lj 


Introducing W(s^) in (7) yields: 
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i = eC 


00 

,1 


£ . -£ 

exp C - \ T F) +1 


-1 mV 


2lV 


r - „ -At Cosh2TT (A+y) -Cosh27r (A-y) 
3 l e l i o J * Cosh2irfX+u)+2Tra 


exp 


(e.-e ) 
1 o 

kT 


d£. 

i 


( 10 ) 


The activation overpotential, p, is included in e^. In a numerical 
solution, made by Bockris and Matthews 5 , the barrier parameters were 
taken as those which produced a consistent Tafel line and was also 
consistent with the variation of the H/T separation factor with 
potential 5 . 

2. The Electrostatic Fluctuation Approach 10 

Using the Fermi Golden Rule of second-order perturbation theory, 
a general expression for the transition probability may be written as: 

W if (£ i } = r- Av Z f|'t'|V'F i dv| 2 5(E. - E f ) , (11) 

f • 

where and ^ are the initial and final state wave functions, v is 
the perturbation operator, and 6(E^ - E^) is the dirac-delta function. 

For any electron transfer reactions, in the adiabatic approximation, 
the initial and final state wave functions can be written as the product 
of the electronic and the solvent wave function in each of these states 
respectively. Thus: 

T =«J>(r,R)x(R) > (12) 

where <f> is the wave function of the fast sub-system, i.e., the electron, 
and x is wave function of the slow sub-system, i.e., the solvent 
molecules, r and R are the electronic and solvent co-ordinates, respec- 
tively. 

Thus, equation (11) can be written as: 
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iiT 2 tt « 
w if = h - Av 


(r , R) x* (R) V (r , R) Xi (R) (Jk (r , R) drdR | 



E^. (13) 


Since the wave function of the heavy particles is localized in space, 
whereas the wave function <{)(r,R) and the perturbation energy 
v(r,R) are distributed over the whole volume, we can write equation 
(13) as: 


"if ' r Av £ljW r l 


Xf X i dRr6(E f 


E i) 


(14) 


Defining L = 


< J ) £ V< I ) jdr | , we get: 


- IT 1 2/ " $ 


X|X i dR|6(E £ - E t ) 


(15) 


Taking the Fourier expansion of the delta function, equation (15) may be 
rewritten as: 


"if‘ e i> - h I I2 

h 


l Av j : 

-A • 


Xf|exp j^(H f - H i )t 


dt. 


(16) 


To solve this matrix element, we have to define H- and H^., which are 

if 

Hamiltonians of the initial and final states. Thus, we need a Hamil- 
tonian for the initial and final states. The Hamiltonian may be 
written as: 


H = H + 
a e 


H + 
s 


V 


e,s 


where a = i,f, H g is the Hamiltonian for the kinetic energy of the 

electron, H is the Hamiltonian for the polar solvent and V is the 
s e,s 

interaction between the electron and the solvent. 

Levich 9 has suggested the Hamiltonian for the solvent H g may be 
written as: 
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The next problem is then to evaluate V . The electrostatic 
approach 10 considers the hydrated ion to be surrounded by dipoles. 

We may divide the dipoles of the solvent into two groups: (1) molecules 
which are arranged closely around the ion and form the so-called inner 
hydration sphere; and (2) all the remaining solvent continuum. However, 
the electrostatic fluctuation approach neglects the interaction between 
the ion and the inner solvation sheath, because the energy of inter- 
action for the ion and the inner layer is such that for these bonds 
hun » kT. Thus, this approach, instead of looking for a Hamiltonian 
comprising of specific interactions, with the ion and the inner sphere, 
considers only the long range interaction between the electron in 
the ion and the continuum solvent. This interaction between the 
electron and the far off dipoles in the solvent is small. However, 
from time to time, a fluctuation from the normal distribution will occur 
(see Section IV) . It has been assumed that such fluctuations occur 
with a frequency sufficient to rationalize observed rates. Using 
macroscopic dielectric theory, we can therefore write: 



D (P. + P . ), dV 

e v in non-m J long 


) 


(18) 


where D is the induction of the electron field, P. ®d P . represents 
e m non-xn r 

the inertial and non-inertial part of the polarization, the subscript 
"long" indicates that the interaction is only with the longitudinal 


JO* 




component of the wave. 
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Now, solving (16) using a harmonic oscillator wave function 11 and 
averaging the matrix element using Feynman's method, one obtains: 


W if ^ - M : 


77 

X, 

"2 

eX P “ 

f ( J - - J . + E ) 2 ) 
f 1 s 


4E kT 

s 

J 


s 


(19) 


where and are the ground state energies of the initial and final 
states and E s the repolarization energy is defined as: 

E s = ( h<0 o /2) Z «,'f - ,W) 2 , (20) 


where u) q is the librational frequency of the solvent dipoles . The 
quantity E g is a crucial quantity of the electrostatic continuum theory. 
However, it cannot be evaluated from equation (20) which is the quantum 
mechanical solution. Thus, in EF model, the basic quantity E g is evaluated 
from classical macroscopic dielectric theory. Thus: 


E = 
s 


hco 


[<& - <&} 

2 2 
= (ne) Z 

f 

1 1 

r-H 

I-H 



e ^ e ^ 

2a, 2sl~ R 

U J 

k 



opt stat 

{ 1 2 J 


( 21 ) 


where a^ and a^ are the radius of the ions between which electron transfer 
is occurring and R is the distance of closest approach. The current can 
therefore be written as: 


i = 2eC 


ti 1 


W if (e i )n(e i )p(e i )de i , 


( 22 ) 


where C is the surface concentration of reactants, n(e.) is the Fermi 
H 1 

distribution and p(e^) is the density of states. Introducing the 

expression for W^f( e ^) (22), we get: 

K1 f exp{-(e - e° ,) 2 /4E kTj 

1 ■ 2 V*Hr 6, - n,) J ' g 5 g n q - -"^yHr rir d i - t23) 
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where : 


O _ ,CK 

e , = e v + E + (J~ - J.) - eri, 
nn' F s k £ x 


where e p is the Fermi energy, and rj is the overpotential; k is the 
transmission coefficient. 

Under the special condition | eri - (J pn , - J?^) | <<£2 > when equation 
(23) is solved by Laplace's method (i.e., by the steepest descent method 
for the function of a real variable), it yields 9 the equation: 


f 


i = eC p*oo K(n,n')kT exp 
H 0 


CJ4 


4kT 


2kT 


J?) 

1 + eri 


2kT 


(24) 


which is of the same form as Tafel's equation. 

Thus, the important point is that the well known Tafel relation 
so widely found to be applicable is obtained from the EF model only when: 


fen - (J f - J?)|«E S . 


This is a very restrictive condition . In fact, the region of over- 
potential where the above condition holds is about 0.4v., whereas, 
experimentally, the linearity continues over a range of >lv. 


IV. COMPARISON WITH EXPERIMENT 

1. Relation between k, . and k, 

het horn 

An argument which has been used to support the continuum electro- 
static theory 12 is that the rate constant for the heterogeneous reaction 
can be predicted from the rate constant for the homogeneous reaction. 

The prediction is: 
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fk. 

hom 


pwl 

z hom 


z het 


where k^ om and \ e t are t * ie rate constants for the exchange reaction 

in solution ( e.g Fe^ + Fe -* Fe + Fe^ ) and at an electrode 
3+ 2+ 

(e.g., Fe + e(M) -* Fe ). However, it is simple to demonstrate 1 3 
that such a relation is expected independently of model, so long as 
the mechanism of both homogeneous and heterogeneous reactions is the 
same. 


2 . Predictions of the Free Energies of Activation from a Continuum 
Electrostatic Theory 

From (23), the free energy of activation is: 


AF* _ 
calc 


(J, 


J.) + E 
x J s 


4E kT 
s 


E s being defined as in equation (21) . 

AF*ai c va lues were computed for 52 reactions, including both the 

2 + 3+ *3+ 2 + 

electron exchange (e.g.,, Fe* + Fe -> Fe + Fe ) and usual electron 
2+ 4+ 3+ 3+ 

transfer (e.g., Fe + C e Fe + C e ) reactions. They were made 
for various ligands, as well. 

The A F * alc is plotted against the experimental AF* value in Fig. 1. 
There is no co-relation between the prediction of the electrostatic 
fluctuation approach with experiment. 



For ions of radius about 3&, water molecules no longer orient 
around the ion, or remain attached to it for significant times. The 


continuum viewpoint (which discounts the effects of such an oriented 

AF°* 

C3.X C 

layer) might then be more applicable. Hence, a plot of - — j — would 

✓ AF° . 

expt 

be expected to tend to unity with increasing radius. The test (Fig. 
II) does not support the applicability of EF. 


3. Variation of 3 with overpotential and the Tafel relationship from 
both approaches 

In Figure III the transfer coefficient, calculated according to 
the electrostatic fluctuation theory, is seen to vary linearly with 
overpotential. Using these values of the transfer coefficients for 
E s = 2ev, the predicted Tafel line is shown in Fig. III. There is no 
linear region, in marked variance with the predictions of the models in 
electrode kinetics. 

In the collisional activation model 24 3 is: 



1 1 - {(<AH°*> - (1 - 3)Fn)/D 1 ) 1 ' 2 
2aX o 1 - {(<AH°*> - A q + Fri)} 1 ' 1 


(26) 


Here, 3 is at first constant, but varies at high overpotential. The 
variation is shown in Fig. III. The agreement with experiment is 


reasonable. 
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V. DISCUSSION 

The fact that the current -potential relation obtained from the 
EF model arises from the parabolic nature of the potential-energy 
profile assumed in that model. Sufficiently small perturbations in 
libratory oscillations of the solvent molecules, far away from the ion, 
must be harmonic, so that introduction of anharmonicity to the present 
model is not acceptable. 

Polaron theory in solids, where a system of interacting oscillators 
are coupled to one another, has postulated fluctuations of the type 
discussed. In a liquid without appreciable periodicity, the applicability 
of the polaron theory is more difficult to accept 14 . Damping of the 
wave would have to be introduced and there is no way at present of 
discovering these factors. 

The rotation-vihrational energy levels in H^O* ion a re not far 
apart compared with the observational significance in electrode kinetics 
(lmev) . As a consequence of the frequent collisions between the solva- 
tion ion and the surrounding water molecules, they are smoothed out. 

The continuous relation observed in electrochemical kinetics between 
rate and overpotential could not be obtained if they were. The opinion 
that thermal activation would produce a structured Tafel line was a 
principal reason for suggesting a model basically different from collision- 
al activation. According to Falk and Giguere 15 , for liquid 1^0 there is 
a continuum of energy states from 60 to 3444 cm ^ (cf. O’Ferrall et al 16 ). 
Thus, classical (continuous) modes of heat transfer are available in 
water and, in particular, aq. H^O*. The appreciable partial molar heat 
capacity (C^) at room temperature for H^O*, indicates a range of frequencies 
giving classical contributions in this molecule (and hence in H^0 + ^ an d 
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provides compelling evidence for the above conclusions. In the gas 
phase, translational energy is converted into vibrational energy 
(after collision) , and it is reasonable to assume that the same 
process occurs in liquid water. The translational energy of free 
waters obeys a Boltzman distribution of energy and is in equilibrium 
with the H^0 + ion in solution. The vibration-rotational energy levels 
of the H^0 + ion will also have a Boltzman distribution (at least up 
to a certain energy level. There are 17 sufficient 0-H bonds in an 
activated state to sustain the reaction over the range of rates observed 
for electrode reactions on various metals. 

If the linear spheres of the reactant and product are different, 
there must be inner sphere activation to satisfy the Franck-Condon 
restriction of the electron transfer process 18 . Thus, the reactant 
and product activated state must be in resonance and identical. 

If the inner sphere is different, activation must occur in it before 
the right activated state can be obtained. The explicit omission of 
considerations of the inner sphere 9 in EF implies that the inner sphere 
for reactant and product have the same energy, without activation. 
Applied to HgO + + H ^ s , such a concept is unreasonable and a model 

based on it 9 unacceptable for such reactions. 

The Born-Landau Hamiltonian arises when the electronic adiabatic 
approximation is used to treat the problem of an electron trapped in a 
polar solvent. This approximation can be applied if the binding energy 

of the trapped electron is smaller than the binding energy of the medium 
electrons. It is a satisfactory approximation in a polar crystal where 
the binding energy of the trapped electron has been shown (Markham and 
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Seitz 19 and Frohlich 20 ) to be of the order of 0.1 ev, much less than 

the binding energy of the medium electrons. With polar liquids , 

however, the situation is different. The binding energy of the 

medium electrons (4-5 ev) is of the same order of magnitude as the 

binding energy of the trapped electrons (1-2 ev) ; the electronic 

adiabatic approximation thus becomes invalid. To calculate the 

binding energy of electrons in polar liquids, Jortner 21 suggests that 

the trapped electron and the medium electron should be treated on an 

equal basis and uses the independent particle treatment, treated in 

terms of self-consistent field theory. The difference between the 

results of this approach and those obtained with the electronic 

adiabatic approximation, is due to the fact that in the self-consistent 

field scheme, the electronic polarisation does contribute to the binding 

energy of the electrons. Thus, the potential becomes of the form 
2r , \ _ 2 


e 

2r 


1 - 1 - 

e 

s 


instead of the 


2r 


,e e „ t 

( op stat 


involved in the expression 


derived from the electronic adiabatic approximation. The reorganisa- 
tional energy would no longer be given in terms of the Born-Landau 
equation. It would correspondingly, not be given by the time-average 
Born energy, but by an equation for the change in the orientation 
energy on ejection of a charge into a polar dielectric. 

Hitherto, we have discussed the different degree of agreement with 
experiment achieved between the two models, and also some fundamental 
difficulties both in the origin of the f luctuational model aid its use 
of the Born-Landau equation. However, we have not discussed the most 
important difference between the two approaches. This is the nature of 
the activation process. The thermal view considers that vibrational 
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excitation, due to collisions between the ion and solvent very near 
to the electrode creates the activated state necessary for electron 
transfer. It does not, thus, differ from well-known concepts in 
reaction kinetics. The electrostatic fluctuation assumes that an 
electrostatic fluctuation in the ion-solvent interaction is the 
origin of the activation. Let us calculate which mode of activation 
gives the higher probability of obtaining activated states of certain 
values. 

Consider a situation where an ion (including the diameter of a 
water molecule in the inner solvation sheath) is r. and a charge of 

ft 

+e is surrounded by N water dipoles. For non-interacting dipoles, 
the probability of the dipole having an angle between 0 and 0 + d0 at 
a distance between R and R + dR from the ion, is: 


P. = 


N 


3 

- exp 

epcos0 

9 

R 2 dR.sin0d0 

2r 

£ R^kT 


CO 

s 



IT 


f 00 3 


o r. 


,, 3 

2r u 


exp 


eycos0 


(27) 


6 R 2 kT 
s 


R dR.sin0d0 


Equation (27) is obtained for non-interacting dipoles. There are of 
course very strong interactions in a polar liquid such as water. 
Following Kirkwood 22 , let us represent the case in water also as a 
system of non-interacting dipoles having the effective dipole moment 

y eff given as: 

y eff = y + g (Ucosy) 

= y(l + g cosy) , 
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where g is the number of neighbour water molecules around any water 
molecule, and cosy is the average cosines of the angles between dipole 
moments of the central molecule and those of its neighbours. In 
evaluating gcosy, Kirkwood 22 considered the nearest neighbours and 
Pople 23 extended Kirkwood's 22 treatment and considered the contribution 
of both first and second layers of water molecules and obtained 
P e ££ = 1.53 as the effective dipole moment of water in water. 

Therefore, in liquid water we consider a system of N-non interacting 
dipoles, however with an effective moment different from that in the 
vapour phase. 


On this basis, eqn. (27) becomes: 


3 

ey eff cos0 

R 2 dRsinQd6 

x ex P 
2r 

w 

e R 2 kT 
1 s 

TT T _ 

00 3 

ey eff cosG 

R 2 dRsin0d0 

x ex P 

K 

o r . 

e R 2 kT . 
s 


(28) 


The average value of the interaction energy between the ion and 
the dipole can then be written as: 


[E] = 


17 L 

r 3 

ey eff cos0 

? 

R dRexp 

* * 
e y effCos0 

J J 2r 3 

or. a) 
i 

e R 2 

c 

e R 2 kT 

c 





sin0d9 


(29) 


7T r 

fCO 


exp 


2r 

o r^ a) 


e^fCos 0 


e R 2 kT 
s 


irdRsined 0 


It has been shown that 1 the integrals can be evaluated to give: 

2 .2 




1 e y eff 


N 2 3 

e r r.kT 
s w 1 


( 30 ) 
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The total interaction energy is then: 


[u] = N [E] = - 


2 2 

6 y eff 


2 3 

e rr.kT 

S 03 1 


( 31 ) 


Following a similar procedure, we can evaluate the average square of 
the energy. This is:* 


2 * 2 4f 

1“ ] ■ tHF 

err. 

S 03 1 


(32) 


The mean square deviation then is given as: 


2 r 2 1 r ,2 

a = [u ] - [u] . 


(33) 


Numerical evaluation of a shows 17 that it is very much smaller than 
2 

[u] . Thus, we can use the Gaussian distribution to evaluate the 
possibility of fluctuation, i.e.. 


H 3° 


exp 

(27 0 2 ff 




2cr 


(34) 


For the proton-discharge reaction, a typical measured activation energy 

at the reversible potential is 20kcals/mole. Thus, we need a fluctuation 

-12 

of (u - [u]) = 20kcals/mole = 1.39 x 10 “ ergs/molecule. Thus, using 

-41 

equation (34) , P comes out to be - 10 
H 3° 

The corresponding probability of the thermal activation can be 
obtained f rom. the Bolt zman distribution and for the 0-H* bond and the 


* The value of the dielectric constant to be used in these expressions 
is that for a time average evaluation, i.e., 80. 
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activation energy of 20kcals/mole. It is 10 17 . Hence, the CA 

approach gives a far higher probability of creating an electric state 
suitable for electron transfer at an electrode than does the 
continuum electrostatic approach. 

So far, we have estimated the probabilities at the reversible 
potential. However, it is important to estimate how these probabilities 
vary as the potential changes. For the hydrogen evolution reaction, 
the activation energy at any overpotential is given as: 

AH^ = (AtfS -3r,F. 

n = 0 

Thus, the variation of the activation energy with overpotential 
is known. We can estimate the change of probability with potential 
using a Boltzman distribution and the f luctuational probability using 
a Gaussian distribution. The results are tabulated in Table I. 

The ratio of the two probabilities are plotted in Fig. V. The conclusion 
is very clear that in most of the experimentally accessible overpotential, 
the thermal probability is much larger than the f luctuational one. 

Thus, the approach arose as a reaction to some misunderstanding 
of the distance between levels in the vibration-rotation spectra of H^O 
in solution. The acceptor states in H^0 + arising from the normal Boltz- 
man distribution outweighs the f luctuational probability over much of the 
range of experimentally obtainable overpotentials, at least for hydrogen 


evolution. 
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TABLE I 


n in 
volts 

AH* 

in kcals/mole 

Probability 

Thermal 

Probability 

Fluctuational 

Probability 
ratio (T/F) 

log 

(T/F) 

0 

20 

O 

1 

4^ 

io" 42 

io 28 

28 

0.2 

17.7 

io- 12 - 4 

10-32. 8 

io 20 

20 

0.4 

15.4 

10 " 11 

IO" 25 

10 14 

14 

0.6 

13.1 

io" 9 * 5 

IO' 17,8 

io 8,3 

8.3 

0.8 

10.8 

10" 8 - 1 

10-12.! 

IO 4 - 3 

4.3 

1.0 

8.5 

io" 6 * 7 

io" 7 * 6 

10>- 4 

1.4 

1.2 

6.2 

io- 4 - 5 

io" 4 - 2 

10 0 - 2 

0.2 

1.4 

3.9 

IO" 3 - 0 

io 1 * 68 

io- 0 - 9 

-0.9 

1.6 

1.4 

IO" 0 - 9 

io" 0 * 4 

io" 0 - 5 

-0.5 



v'-v 










cole 
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ANNOTATIONS FOR FIGURES 


Figure I. The free energy of activation for redox reactions in 


solution. 


AF 


represents the value obtained by 

exptl 


equating the experimental rate of Z exp 


AF 0 *" 1 


RT 


is the calculated bimolecular collision number; AF 
is the value obtained from the electrostatic solv 

fluctuation viewpoint. There is no correlation. 


, where Z 
o* 


. AF 


Figure II. 


theoret 


as a function of the sum of the radii of the 


AF 


exptl 

ions. The value does not tend to unity at high radii. 


Figure III. S as a function of potential , electrostatic and thermal. 

That calculated on the electrostatic model with a reorganisa- 
tion energy of 2ev varies much more over the range O-lv 
overpotential than does that calculated from the electro- 
static view. Experimentally, the coefficient is nearly 
constant with potential at over at least 0.5 volts. 


Figure IV. The experimental Tafel line for overpotential as a function 
of log i in the system: H^0 + + e H 2 [Points ) . Broken 
lines: thermal model. Continuous line: electrostatic 
fluctuation theory. 
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Figure V. The relative probability of the activation of protons into 
states above the ground state by amounts varying from lev 
to 0. P T is the probability arising from the theory in which 
there is equilibrium between the heat sink and the surrounding 
solvent. Pp is the probability calculated on a fluctuation 
of the electrostatic energy of the ion. 
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MAJOR POINTS FROM SEN'S WORK 

1. The relation of friction to potential when the two interfaces 

are metal-solution, can be interpreted in terms of the 
repulsion of two Gouy layers. 

2. The BDM isotherm allows predictions of a capacitance minimum. It 

is not inconsistent with thermodynamic reasoning. 

The polarizability term in the double layer equations may 
become important at high enough fields. 

3. The polaron theory of electrode processes is inconsistent with 

line broadening in solution. 

4. The polaron theory of electrode processes is not consistent with 

the trend of electrochemical facts about proton transfer. 

5. The polaron theory is not able to give the Tafel line. 

6. The limits of the Gamow approximation for the probability of trans- 

fer are given. 

7. The Gurnean and Weissian approaches to electrode kinetic 

formulations are given. 

8. BEBO gives good results for proton transfer. 






